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ABSTRACT 


This article presents the analysis and design methodology of a steel frame 
with eccentric diagonals in the shape of an inverted V (Chevron) with the link 
element in the center. The analysis and design procedure is illustrated using the 
CEINCI-LAB Computer System and a detailed explanation of the new functions is 
illustrated. Through these functions the user can perform seismic and static 
analysis step by step, with this the calculation procedure can be understood. 


An example of a pre-design of a steel frame with eccentric diagonals is 
illustrated, the seismic analysis is carried out to obtain displacements, drifts and 
floor shear. In addition, the stability of the structure against overturning is verified. 
A static analysis is carried out with the states and load combinations of the 
Ecuadorian Construction Standard 2015 (NEC-15), this allows defining the demand 
“D” in the elements. The capacity "C" of elements is calculated and the relationship 
"D / C" is presented graphically so that the user can verify that this factor is less 
than 1. Finally, the analysis by capacity of the structure is carried out to define the 
design. 


keywords: CEINCI-LAB Computer System. Seismic Analysis. Static Analysis. 
Steel structures. Eccentrically Braced Frames. NEC - 15. Steel Design. 


NUEVAS FUNCIONES DE CEINCI-LAB PARA EL 
ANALISIS Y DISENO DE PORTICOS DE ACERO CON 
ARRIOSTRAMIENTOS EXCENTRICOS 


RESUMEN 


En este artículo se presenta la metodología de análisis y diseño de un 
pórtico de acero con diagonales excéntricas en forma de V invertida (Chevrón), 
con el elemento de enlace en el centro. Se ilustra el procedimiento de análisis y 
diseño mediante el Sistema de Computación CEINCI-LAB y se describe una 
explicación detallada de las nuevas funciones. Mediante estas funciones el usuario 
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puede realizar análisis sísmico y estático paso a paso, con ello, se puede entender 
el procedimiento de cálculo. 


Se ilustra un ejemplo de prediseño de un pórtico de acero con diagonales 
excéntricas, en el cual, se realiza el análisis sísmico para obtener 
desplazamientos, derivas y cortante de piso. Además, se verifica la estabilidad de 
la estructura ante volteo. Se realiza un análisis estático con los estados y 
combinaciones de carga de la Norma Ecuatoriana de Construcción 2015 (NEC- 
15), esto permite definir la demanda “D” en los elementos. Se calcula la capacidad 
“C” de elementos y se presenta de manera gráfica la relación “D/C” para que el 
usuario pueda verificar que este factor sea menor a 1. Finalmente, se realiza el 
análisis por capacidad de la estructura para definir el diseño. 


Palabras Claves: Sistema de computación CEINCI-LAB. Análisis Sísmico. 
Análisis Estático. Estructuras de Acero. Pórticos Arriostrado 
Excéntricamente. NEC — 15. Diseño de Acero. 


1. INTRODUCCIÓN 


Los primeros estudios de pórticos arriostrados excéntricamente se 
desarrollan en Japón y en Estados Unidos. Trabajos como los de Fujimoto et al. 
1972; Tanabashi et al. 1974; y Roeder & Popov, 1978 sentaron las bases para el 
estudio y el uso de estos elementos, y es así como, en 1980 se desarrollan 
numerosos estudios sobre el comportamiento no lineal del elemento de enlace. 
Además, se estudia el comportamiento cíclico de estos elementos en trabajos 
como los de Manheim & Popov, 1983; Hjelmstad & Popov, 1984; Malley & Popov, 
1984; Kasai & Popov, 1986; Ricles & Popov, 1989; Engelhardt & Popov, 1992. 


Estudios experimentales como los llevados a cabo por Yang, 1985; 
Roeder et al.1987; Whittaker et al., 1989, permitieron que se incluya el uso de 
diagonales excéntricas en disposiciones como la UBC 1988 y posteriormente en 
las disposiciones sísmicas del AISC. 


A lo largo del tiempo se han realizado diversas investigaciones sobre 
estos elementos, entre ellas, se destacan trabajos como Mansour et al., 2011, 
Berman 8 Bruneau, 2008, dónde el enfoque se basa en la diversidad de formas 
que puede tener un elemento de enlace, incluyendo en estos, elementos con 
forma I, secciones construidas en base a placas y secciones tipo caja (cuadrados 
o rectangulares). Posterior a ello y con el cambio de filosofía en el diseño de 
estructuras sismorresistentes, de un diseño basado en capacidad, a un diseño 
basado en desempeño, surge la idealización de enlaces reemplazables como lo 
exponen autores como Ramadan & Ghobarah, 1995; Mansour, et al. 2011; 
Dusicka & Lewis, 2010. 


En Ecuador es común emplear pórticos duales como sistemas resistentes 
a cargas laterales; los pórticos duales son compuestos por pórticos resistentes a 
momentos y en algunos ejes se incluye elementos diagonales, que pueden ser 
concéntricas o excéntricas en función de los criterios del arquitecto e ingeniero 
responsables del proyecto. En ese artículo se estudia el análisis y diseño de 
pórticos duales con arriostramientos excéntricos mediante CEINCI-LAB. 
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2. COMPORTAMIENTO DE PÓRTICOS ARRIOSTRADOS 
EXCENTRICAMENTE 


Las diagonales excéntricas son aquellos elementos que no logran unirse 
entre sí, sino, que permiten un espacio entre ellas, a este espacio se lo conoce 
como excentricidad (Aguilar, 2015). 


El diseño de estructuras con arriostramientos excéntricos se utiliza 
principalmente en edificios de alta esbeltez (relación de la altura de un edificio 
debe mucho mayor que 4 veces su ancho), en los cuales se busca disipar energía 
a través de la fluencia del elemento de enlace. La ubicación de dichas 
excentricidades está ligada tanto con la ductilidad, como la rigidez elástica de la 
estructura, por tanto, se debe evitar colocar estos elementos de enlace cerca de 
columnas (Figura 1a), puesto que, se pueden generar rótulas plásticas en la parte 
superior de las columnas y producir una falla de la estructura por piso flexible 
(Aguilar, 2015). 


La tipología estructural de pórticos arriostrados excéntricamente se 
compone de columnas, vigas y riostra diagonal. Además, existe una característica 
distintiva en estos pórticos y es que, al menos uno de los extremos de cada riostra 
diagonal está conectado de manera que la fuerza de la riostra diagonal se 
transmite a través del corte y la flexión del elemento enlace (link element). En la 
Figura 1 se indican varios tipos de pórticos con excentricidad horizontal entre los 
puntos de intersección de las dos líneas centrales de las riostras diagonales con la 
línea central de la viga. 


A Je 

















(a) (b) (c) 
Figura 1 Pórticos con diagonales excéntricas (basado en Popov et al. 1987). 


Los pórticos con arriostramientos excéntricos pueden proporcionar una 
rigidez elástica comparable a los pórticos arriostrados concéntricamente, sobre 
todo cuando se utilizan enlaces cortos (Roeder y Popov, 1978; Libby, 1981; 
Merovich et al., 1982; Hjelmstad y Popov, 1983; Malley y Popov, 1984; Kasai y 
Popov, 1986a, 1986b; Ricles y Popov, 1987a, 1987b; Engelhardt y Popov, 1989a, 
1989b; Popov et al., 1989). Además, puede generar ductilidad y capacidad de 
disipación de energía en el rango inelastico; siendo comparable a la de los 
pórticos especiales resistentes a momento; siempre que los enlaces no sean 
demasiado cortos (Comentarios AISC 341). 
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En el trabajo de Hjelmstad & Popov (1984), se estudia la variación de la 
rigidez lateral de un pórtico con arriostramiento excéntrico, respecto a la rigidez del 
pórtico sin diagonales, estas variables varían en función de valores de la longitud 
del enlace y para distintos valores de longitud de vano. El pórtico de la Figura 2a 
es el que demuestra a través de estas variaciones, mostradas en la Figura 2b, que 
existe mayor rigidez lateral cuando se emplean enlaces cortos, lo cual, resulta 
ventajoso para el control de la deriva de piso. La variable Kpy, es la rigidez lateral 
del pórtico arriostrado excéntricamente y la variable Kpry se refiere a la rigidez 
lateral del pórtico resistente a momento. 


k He~ 
p | 
h 
l 
¡A Í 
(a) 


Porticos Arriostrados Excéntricamente 
Enlace en el vano central 


Rigidez relativa K. 


eL 
(b) 

Figura 2. Resultados del pórtico analizado: a) Pórtico de estudio; b) Variación de 
la rigidez lateral entre el pórtico con arriostramiento excéntrico versus la 
rigidez lateral el pórtico a momento. 

Fuente: Este estudio. 


Acorde a Popov et al. (1987) existe una relación entre la fuerza cortante en 
el enlace V;;,, con el cortante estático de diseño de cada piso, acumulada desde 
la parte superior hasta el nivel correspondiente de una estructura Vp;,,. La relación 
mencionada se puede expresar mediante la ecuación 1, la cual se puede obtener 
en base a una aproximación del equilibrio de fuerzas en el extremo A del pórtico 
de la Figura 3. 
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Figura 3. Aproximación de diagrama de cuerpo libre (basado en Popov et al. 


1987). 


La ecuación 1 es una aproximación para obtener la fuerza en el enlace a 
través de variables conocidas. Mediante la modelación de una gran cantidad de 
pórticos con diagonales excéntricas, se puede visualizar la relación entre la fuerza 
cortante en el elemento de enlace Veyrace y el cortante de diseño en cada piso 
Vprso, en función de la geometría del pórtico (h/L), y la relación entre la longitud 
del enlace y del vano del pórtico (e/L). Estas relaciones se expresan en la Figura 
4, donde se puede notar que la relación entre los cortantes es mayor cuando la 
relación e/L es menor, además existe una relación proporcional con respecto a 
h/L. Además, se plantea la ecuación 2 mediante regresiones multivariables que 
describe las curvas de la Figura 4, la ecuación 1 es conservadora y no depende de 
la relación entre e/L que es influyente en la rigidez de la estructura. En la Figura 4 
las líneas continuas son los valores determinados al modelar los pórticos y las 
líneas punteadas representan las calculas con la ecuación 2, mientras que la línea 
de color negro representa la calculada con la ecuación 1. 


Venlace'Vpso 
= = o 
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al =00% 
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Figura 4. Variación del cortante en función de la relación de aspecto y longitud del 
enlace. 

Fuente: Este estudio. 
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2.1. Filosofía de diseño en pórticos de acero con 
arriostramiento excéntrico 


La filosofía de diseño se fundamenta en la fluencia del elemento de enlace 
debido a corte, flexión o ambos efectos, en función de la longitud de este; además 
se espera que, las vigas fuera del enlace, las diagonales y las columnas se 
mantengan elásticas, aunque es posible aceptar incursiones inelasticas sin 
comprometer la estabilidad de la estructura (es decir fluencia en estos elementos). 
En el diseño se asume que los enlaces podrán desarrollar su máxima capacidad 
esperada y bajo esas fuerzas se diseñan los otros elementos. En la Figura 11 se 
muestra la deformada de un pórtico con arriostramientos excéntricos. 


Enlace (Link) 

e  Fluencia por flexión 

e Fluencia por corte 

e  Disipación de energía 


Conjunto elástico 

e Viga fuera del enlace 
e  Diagonales 

e Columna 





Figura 5. Deformada de un pórtico ante cargas laterales. 


Para el análisis estructural, en el diseño por capacidad, es necesario 
reemplazar el elemento de enlace por las acciones máximas esperadas debido a 
la fluencia de estos elementos (V yp hace referencia al cortante máximo probable y 
M mp se refiere al momento maximo probable); este estado de cargas sustituirá a 
la combinación por sismo en esta etapa de análisis. En la Figura 6 se indica el 
equilibrio de fuerzas debido a la fluencia del elemento de enlace. 











Figura 6. Fuerzas debido a la fluencia del elemento de enlace. 
El tipo de falla del enlace depende de la longitud de este, respecto a la 
longitud critica. La longitud critica o excentricidad critica es la longitud del enlace 
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que permite desarrollar el corte y momento plástico en la sección de manera 
simultánea (V, y M, respectivamente). Esto se demuestra al realizar el equilibrio 
de fuerzas del diagrama de cuerpo libre de la Figura 7a, tal y como se ve en la 
Figura 7b. 





CIEI 
V 
M 
(a) 
j e < 
M, D L 
Ma V, 


(b) 
Figura 7. Fuerzas y momentos en elementos de enlace; a) Diagrama de cuerpo 
libe y diagrama de corte y momento; b) Equilibrio plástico. 


Al comparar la longitud del enlace respecto al parámetro de longitud critica 
se puede estimar que la fluencia a corte se presentará en elementos cortos. Es 
decir, longitud efectiva menor que la longitud crítica y en caso contrario la fluencia 
será debido a flexión en el enlace, esto se ilustra en las Figuras 8a y 8b. 


La fuencia por corte ocurre 


‘ La fluencia por fexidn ocurre 
v ¿ u cuando 


u | n cuando 
M-<M, y V=YW,, = l MM, y V<V, 
= 2M, 7 2M, 
Yy V 


(a) (b) 
Figura 8. Comportamiento de elementos de enlace en función de la longitud; a) 
Fluencia por corte; b) Fluencia por Flexión. 
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El comportamiento del elemento de enlace, expresado en términos de 
rotación inelástica máxima y capacidad a corte (Figura 9), depende del cortante y 
momento plástico de la sección de este elemento, en base a la relación expresada 
en la ecuación 3. 





e 
Factor = ——— (3) 
Mp/Vp 
I ' 
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'Fuencia, | flexion 
d rn Y 
à Corte+ j 
ı Flexión ı 
4 t 
1.6 26 
(a) 
Vn = Vp 
Vn 2Mp 
€ 


(b) 
Figura 9 Comportamiento del elemento de enlace; a) Rotación inelástica; b) 
Relación de cortante. 


3. DISEÑO CON DIAGONALES EXCÉNTRICAS 


Para una mejor comprensión de la metodología de diseño de elementos de 
acero, es necesario conocer las ecuaciones básicas para determinar la capacidad 
axial, a corte y flexión de los elementos que componen un pórtico de acero. 


Los artículos de Cagua et al. (2021) y Cagua et al. (2021,2) presentan la 
metodología de prediseño, análisis y diseño de pórticos a momento y de pórticos 
con arriostramiento respectivamente. Esta sección presenta un resumen de los 
principales aspectos de estas temáticas. 


3.1 Capacidad de carga axial y flexión 


La capacidad axial (P,,) de un elemento viene dada por la expresión de la 
ecuación 4, donde Fcr corresponde al esfuerzo crítico de Euler y Ag representa el 
área transversal de la sección. La capacidad axial está relacionada con la relación 
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de esbeltez del elemento A, la cual se calcula mediante la ecuación 5, que a su 
vez depende del factor de longitud efectiva K, la longitud L y radio de giro r. Es es 
el módulo de elasticidad del material y Fy es el esfuerzo de fluencia. Cabe 
destacar la importancia del pandeo como factor de generación de inestabilidad 
elástica presente en elementos a compresión y crítica en elementos esbeltos. 


P, = Fo * Ag (4) 
K-L 
Ja (5) 
r 
1? + Es (6) 
FCTguLER = K-I? 
T) 


FCcTfgyrer < 0.44 : Fy > Fer = 0.877 * FCrEyLER 
(7) 


Fy 
FCYeuer = 0.44 : Fy > Fer = 0.658FC"EULER * Fy 


El factor de longitud efectiva K de columnas para pórticos sin 
arriostramientos, se puede calcular mediante la fórmula 8, la cual depende de las 
rigideces rotatorias de los nudos y con sus condiciones de borde, que, para su 
efecto, se calculan con la fórmula 9. 





8 
G4 + Gp + 7.5 (8) 


27), 
xP), 


Donde G, y Gg son las rigideces rotatorias en el nudo inicial y final 
respectivamente, I es la inercia de la sección, L es la longitud del elemento, los 
subíndices c y b se refieren a columnas y vigas respectivamente. Si el extremo de 
la columna tiene un apoyo articulado entonces se puede tomar un valor de G = 10, 
si tiene un empotramiento G = 1 en ese borde. 


pa Pte GAT 


G= (9) 


La capacidad nominal a corte (Vn) se determina con la ecuación 10. 
Vn=0.6:Fy-A,,* Cv, (10) 


Donde Aw es el área del elemento que resiste el corte; en perfiles tipo | se 
asume que Aw es el área del alma y se calcula con la ecuación 11; Cv, es un 
factor debido a la resistencia del corte y se determina con la ecuación 12. 


Ay =(d-2:tf)t, (11) 
= Cuando See 2E 
tw 


Cv, =1 
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= Caso contrario 


(12) 
jije eE 





Co =— > Y7 
a h/tw 


La variable k, esta intimamente relacionada con la presencia de 
rigidizadores transversales en el enlace, por lo que, si se trabaja de la manera mas 
conservadora, como lo hace CEINCI-LAB, su valor es de 5.34, caso contrario 
debe ser calculado con la fórmula 13. 


5 


a 2 
(7) (13) 
ky = 5.34 cuando > > 3.0 


k, =5+ 


Donde a es la distancia libre entre rigidizadores. 


La capacidad a flexión de una sección | con doble simetría y compacta se 
puede determinar con la ecuación 14. 


Mp = Fy: Z; Cuando Lp < Lp 
Lb — Lp 

Lr — Lp 

Fcr : Sx < Mp Cuando Lb > Lr 





Mn =; C: Mp — (Mp - 0.7: Fy Sx)-( )| < mp Cuando Lp<Lb<Lr_ (14) 


Donde C, es el factor de modificación por pandeo lateral-torsional para 
diagramas de momento no uniformes y Fer es un esfuerzo crítico reducido para 
flexión. Las expresiones matemáticas para calcular estos valores se encuentran en 
el capítulo F del AISC 360 -16. 


Las longitudes Lp, Lr se pueden calcular con las ecuaciones 15 y 16. 


|E 
Lp = 1.76: ry: Fy (15) 


E 
0.7- Fy 









Lr = 1.95: 1%. ° 





Sy hg 
2 
Nts = 

Para secciones H con doble simetria y alas rectangulares entonces: 


Iy: hè 
w 7 4 
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Donde ry es el radio de giro de la sección respecto al eje y-y, en este 
caso corresponde al eje débil de las secciones l; J, es la constante torsional de la 
sección; S, es el módulo elástico de la sección respecto al eje x-x; Para secciones 
| con simetría c = 1, h, es la distancia entre los centroides de los patines. 


Los cálculos que se pueden realizar en CEINCI-LAB tienen la 
consideración principal de tomar la capacidad a flexión como aquella lateralmente 
apoyada, donde Lb < Lp; es importante considerar que todas las fibras se pueden 
plastificar, para ello se debe utilizar un acero con comportamiento dúctil. De 
acuerdo con la normativa, los elementos que resistan cargas laterales deben ser 
altamente dúctiles (AISC 341 tabla D1.1) y deben ser arriostradas de acuerdo con 
la ecuación 20. 


Mn = Mp =Fy:Z, (19) 





Lb = 0.095 -ry- (20) 


E 
Ry : Fy 
3.2 Demandas de carga lateral 


Hasta este punto se realiza un resumen de las ecuaciones para el cálculo 
de la capacidad axial, a corte y flexión de los elementos. Adicionalmente, es 
necesario calcular la demanda de los elementos a través de un análisis sísmico 
que permita determinar las cargas laterales y un análisis estático que permita 
determinar las acciones en cada elemento debidas a las cargas gravitacionales de 
la estructura. 


Para la comprobación de la demanda versus la capacidad de los 
elementos, en base a la metodología de Factores de Carga y Resistencia (LRFD 
por sus siglas en inglés, Load and Resistance Factor Design), se debe cumplir 
que: 


O: Ra 2 Y Vi Qi = Qu (21) 

Donde @ son factores de reducción de la capacidad nominal Rn; y; son 

factores que amplifican los efectos de las cargas; Q; corresponde a las cargas 

aplicadas; Q,,, es la carga última o efecto en los elementos con las solicitaciones 
amplificadas. 


La capacidad de elementos sometidos simultáneamente a cargas axial y 
de flexión se determina evaluando la acción simultánea de cargas axiales y de 
flexión, tal como se indica en las ecuaciones 22 y 23, las cuales dependen de los 
valores de Pu y Pn para su selección. 


Para 








Ø: Pn | Ø-Pn 9 O-Mn~ (22) 
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P Pe 269 Wp oO oi (23) 
on > = 
ara @-Pn 2-0-Pn Ø- Mn 





Donde Pu es la carga última aplicada sobre el elemento, Pn es la carga 
nominal que resiste el elemento, Mu es el momento último aplicado al elemento, Ø 
son factores de reducción de la capacidad nominal, y Mn es el momento nominal 
que resiste el elemento. 


En pórticos duales, es decir pórticos resistentes a momentos y pórticos 
con arriostramientos excéntricos, es importante determinar en las uniones de las 
vigas-columnas de los pórticos resistentes a momento que se cumple el principio 
de viga débil y columna fuerte. La ecuación 24 permite el cálculo de este principio 
garantizando que las columnas sean más resistentes que las vigas y 
estableciendo la base del diseño por capacidad. 





T 
> ide < 1.0 
Y Mpc (24) 
Ymi yz 1 (r funy 
c= col - -— 
á z Ag (25) 
Y Mpo =>) (1.1 Ry Fy Zviga + Muyiga) (26) 


Donde Y Mpc es la sumatoria de los momentos plásticos nominales de las 
columnas que llegan a la junta; Y Mpv es la sumatoria de los momentos plásticos 
nominales de las vigas que llegan a la junta. Zcol y Zviga son los módulos 
plásticos de las columnas y vigas. Puc es la carga axial requerida para la columna 
usando las combinaciones de carga para columnas (se debe incluir la sobre 
resistencia); Ag es el área bruta de la sección; Ry, es el factor de esfuerzo de 
fluencia probable (para A36 este valor es 1.3 acorde a la NEC-SE-AC 2015) y 
MU,iga es el momento adicional basado en las combinaciones de carga 
(producido por la fuerza de corte localizado en la zona donde probablemente se 
tendrá la articulación plástica multiplicado por la distancia hasta la junta o nudo). 


4. PREDISENO DE PORTICOS CON ARRIOSTRAMIENTOS 
EXCENTRICOS UTILIZANDO CEINCI-LAB 


En este artículo se realizará el análisis del pórtico correspondiente al eje 1 
de la Figura 10a, en elevación el pórtico es presentado en la Figura 10b, en el 
mismo se incluirá diagonales en V invertida en el vano central, es decir entre los 
ejes B y C; no obstante, los arriostramientos no convergen en un mismo punto, 
sino que dejan una excentricidad de 0.50 m en el tramo central. 


Los pórticos 1 y 2 presentan la misma geometría, por lo tanto, la rigidez 
lateral de cada pórtico es igual y se puede realizar el análisis y diseño de cada 
pórtico de manera independiente asumiendo un área colaborante de 6m para la 
distribución de cargas gravitaciones. 
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Figura 10 Geometría; a) Vista en planta; b) Vista en elevación del pórtico en 
análisis 


En cuanto a las cargas gravitatorias, por facilidad de cálculo el peso propio 
de la losa se asume como 200 kg/m?, dando como resultado que la carga muerta 
para los entrepisos se considera de 500 kg/m? y para cubierta de 350 kg/m?. La 
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carga viva depende de la ocupación, en este caso se asigna 250 kg/m? a los 
entrepiso y 100 kg/m? a cubiertas (tabla 1). 


Tabla 1 Cargas gravitacionales 





























Parámetro Da Observación 
Peso propio de la losa = 200 Losa Deck de 12 cm 
Otros (Instalaciones) = 16 
Enlucido masillado = 44 2 cm de enlucido 
Acabado de piso = 40 Asumido 
Peso propio de elementos 50 Astimido 
estructurales = 
Peso mampostería = 150 Asúmido (Dep ende ge proyecto 

arquitectónico) 

Carga permanente piso = 300 
Carga permanente cubierta = 150 Cálculo anterior sin mampostería 
Carga viva piso = 250 NEC_SE_CG 2015, pág. 25 
Carga viva cubierta = 100 NEC_SE_CG 2015, pág. 25 

















Se asume que la estructura se encuentra emplazada en la ciudad de Quito 
sobre suelo tipo D; en base a ello se estima los parámetros de la tabla 2. 


Tabla 2 Parámetros que intervienen en el cálculo del Corte Basal. 





CARGAS DE SISMO 7 ] 
Parámetro | Variable | Valor Observación | Referencia 

Zona Sismica v - - Tabla 1. Sec.3.1.1 Pág. 27 
Factor de aceleración en la zona Z DA - Tabla 1. Sec.3.1.1 Pag. 27 
Tipo del perfil del suelo D - - Tabla 2. Sec.3.1.1 Pag. 30 
factor de sitio Fa Fa 12 - Tabla 3. Sec.3.2.2 Pág. 31 
factor de sitio Fd Fa 1.19 Tabla 4. Sec.3.2.2 Pág. 31 
factor comportamiento inelás, suelo Fs 1.28 Tabla 5. Sec.3,2,2 Pág. 32 
Factor asociado al periodo de retorno f 1 - Sec 3.3.1 Pág. 34 
Relación de amplificación espectral n 2.48 = Sec.33.1 Pag. Y 
Aceleración en T=To Safa] | 1:19 5,=nZF, Sec.3.3.1 Pág. 34 
Periodo Limite en T=To Tofs] 0.127 do Sec.3.3.1 Pág. 35 
Periodo Limite en T=Tc Tc [s} 0.698 Te =0.55F, Hd Sec. 3.3.1 Pág. 34 
Periodo Limite en T=T. T [s] 2.856 T,=2.4,F, Sec 3.3.1 Pág. 34 
Coeficiente Ct Ct 0.073 - Sec.6.3.3 Pág. 62 
Coeficiente para Calculo de Periodo a 0.75 F $ec.6.3.3 Pág. 62 
Altura total del elemento hn [m] 12.24 . Planos 

Periodo teórico método 1 Tt [s] 0.478 T= CAS Sec.6.3.3 Pág. 62 
Periodo teórico método 1 mayorado | T2[s] | 062 | 12=1.3'T1 | Sec.6.3.3 Pag. 63. 


El espectro para la ciudad de Quito se presenta en la Figura 11, de color 
azul se tiene el espectro elástico de diseño y en color verde el espectro reducido 
mediante el factor R, en este caso se asume un valor de 6 de manera 
conservadora respecto a tener fuerzas laterales más grandes para este sistema 
estructural. La Norma permite emplear valores superiores a 6 para este sistema 
estructural. 
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Espectro de Diseño NEC 15, Amortiguamiento 5% 






—— Especro Elástico 





— Espectro inelaséco R = 6.00 


Aceleración (9) 





Figura 11 Espectro de diseño para Quito — Suelo D. 


La estimación del periodo fundamental de vibración se puede realizar 
mediante la expresión 27, como se recomienda en la NEC-15. 


T=Ct-h2 (27) 


Los coeficientes para pórticos de acero con arriostramientos, concéntricos 
o excéntricos, se utiliza Ct = 0.073 y a = 0.75, con ello se obtiene un periodo de 
0.477 segundos. 


Para este periodo de vibración, la estructura en análisis, tendría una 
aceleración espectral Sa para el primer modo de vibración correspondiente a 
0.198 g, que coincide con la meseta del espectro. El cortante basal se puede 
determinar con la ecuación 28. 


_ EST a 
R : Øp * Oe 

Donde V es el cortante basal; Sa(T) es la aceleración espectral que se 
determina en función del periodo fundamental de vibración de la estructura T como 
se indica en la ecuación 28; R es el factor de modificación de respuesta sísmica; 
Øp y Øe corresponden al factor de irregularidad en planta y elevación, 
respectivamente; W es la carga reactiva al sismo. En la tabla 3 se presenta los 
valores de las variables mencionadas y el cálculo del coeficiente sísmico, el 
cortante basal es 19.84% de la carga reactiva total. 
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Tabla 3 Coeficiente sísmico. 





Parámetro Variable Valor Unidad Obs. Referencia 
Aceleración en T=To Sa, 1,19 g T1 < Tc Sec.3,3.1 Pág. 34 
Factor de irregularidad en planta Op 1 - Tabla 13. Sec. 5.2.3 Pág. 50 
Factor de irregulandad en De 
elevación 1 Tabla 14. Sec. 5.2.3 Pag, 51 
Factor de reducción de R 6 
respuesta - Tabla 16. Sec. 6.3.4 Pág. 65 
Factor de importancia i 1 Tabla 6. Seco. 4.1 Pág. 39 
Coeficiente sismico Coef, 0,1984 - - Sec, 6,3,2 Pág. 61 
PANA k a = Teos Sec. 6.3.5 Pag. 67 


En el caso general, la carga reactiva de sismo se compone del 100% de la 
carga muerta en los pisos que pueden vibrar durante una excitación externa, para 
este caso se asume 100% de carga muerta y un 25% de participación de la carga 
viva, este porcentaje de contribución de carga viva se tomó a consideración de los 
autores puesto que representa en mayor seguridad estructural y se encuentra 
presente en normativas extranjeras como la normativa peruana. Se destaca que la 
norma de Perú E030 (2018) considera un porcentaje de la carga viva como se ha 
realizado en este artículo para el análisis sísmico. En cambio, la NEC-15 no 
considera la carga viva, únicamente lo hace para bodegas y almacenamiento, en 
el cual se toma 25% del valor de la carga viva. 


En la tabla 4 se presenta el cálculo de la carga reactiva para cada nivel, el 
área típica de piso se obtiene en base al área sombreada de la figura 10a. 


Tabla 4 Carga reactiva de sismo W. 


W=CM+ 
% CV [T] 


h [m] CM [T/m?] CV [T/m2] Área [m2] 








Como la carga reactiva es 210.38 T, entonces el cortante basal 
corresponde a 19.84% de ese valor, que es igual a 41.74 T. 


V = %:W = 19.84% : 210.38 29 
V = 41.74T (29) 
En la NEC15 (2015) se indica el procedimiento para determinar las fuerzas 

laterales en cada piso, para ello se aplica la ecuación 30 mostrada a continuación: 


Wj x hj* 


DAA (30) 
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Donde Wj, Wi es la fracción del peso en cada nivel; hj, hi es la altura de 
los niveles en je í, el valor de “k” se relaciona con el periodo de la estructura y 
varía entre 1 y 2, este valor refleja la forma en la que se distribuyen las fuerzas 
laterales en cada nivel de la estructura; para valores de T < 0.5 segundos k = 1, 
para periodos mayores a 0.5 seg este valor se calcula como k = 0.75 + 0.5 *T, 
el valor máximo de k es 2. V es el cortante basal que se distribuirá en cada piso 
como fuerzas laterales Fj. 


Tabla 5 Fuerzas laterales y cortantes en cada piso. 


Vbasal 
— O, 
Vbasal = 19.84% TOTAL = 41.74 [T] 


Se asume k = 1 para la Distribución de Fuerzas en la altura 


H [m] Hi*W [T*m] 1 | F [T] 
206.55 
371.79 


537.03 
468.18 
1583.55 








En la Figura 12 se representa el cortante en cada piso. 


Corte en cada piso (T) 


4 E 12.34 
3 M 26.49 
2 IS 36.29 


ee ee) 41.74 
0.00 10.00 20.00 30.00 40.00 50.00 


Figura 12 Cortante de piso. 


Una vez calculado el cortante basal en base a la normativa es posible 
aplicar la ecuación 31 con las variables del piso 1 y con esto se tiene la relación 
entre la fuerza cortante del elemento de enlace y el corte acumulado del piso. Se 
aplica esta metodología en el prediseño de los elementos, posteriormente se 
determinan la demanda en cada elemento mediante la aplicación del método de la 
rigidez programado en CEINCI-LAB. 
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L ( 5 ns 
cum 
T (o. 075- E — 0.035: E -— 0.05: E) = 0.075) 
pats (-0.5:7 +1) 70 
0.5 0.5 
en .07 s- (=) ‘0.035: (23) - 0.05- (25) - 0.075) (31) 
Viink _ 
Va 0.42 


En la tabla 6 se muestra los valores de cortante para cada elemento 
enlace en cada nivel del pórtico, se asume que el cortante de piso que debe 
resistir los elementos del vano con arriostramientos es 75% del cortante de piso. 


Tabla 6 Cortante en los elementos de enlace de cada nivel. 





E a A 

[T] [T] 

1 3.6 31.3 0.42 13.1 
2 2.9 27.2 0.32 8.7 
3 2.9 19.9 0.32 6.3 
4 2.9 9.3 0.32 2.9 





Para definir las dimensiones de la sección del enlace se emplea la 
combinación de cargas más conservadora para cargas gravitacionales y el efecto 
del sismo (NEC SE CG 2015), esto es: 


Combinación = 1.2D + L + 1.0E (32) 


Hasta este punto se analiza el efecto del sismo, a continuación, se 
requiere analizar el efecto de las cargas gravitacionales para posteriormente 
superponer los efectos. La carga muerta CM en el piso inferior es 500 kg/m? y la 
carga viva CV es de 250 kg/m?. El ancho colaborante perpendicular al pórtico de 
análisis, se asume con un valor de 6 metros (Figura 10a), por lo tanto, la carga 
uniformemente distribuida en la viga con la combinación de cargas se define por: 


q = (1.2: CM + 1.0 : CV) * Leotaborante (33) 


El valor de la carga “q” es 5.1 T/m. El modelo de análisis simplificado para 
la viga con los arriostramientos excéntricos se muestra en la Figura 13a, el 
diagrama de corte y de momento se presentan en la Figura 13b y 13c. 
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Figura 13 Modelo de análisis simplificado para vigas con arriostramientos 
excéntricos; a) Modelo de análisis; b) Diagrama de Cortante; c) 
Diagrama de Momento Flector. 

Fuente: Este estudio. 


Al reemplazar los valores de las variables geométricas y de carga en los 
diagramas de corte y momento se tiene los diagramas de la Figura 14. 


542 Brian Cagua, Roberto Aguiar, Julia Pilatasig y Alejandro Bonilla 





(b) 
Figura 14 Diagramas de corte y momento obtenidos con fórmulas; a) Diagrama de 
Cortante; b) Diagrama de Momento Flector. 


e  Prediseño del elemento de enlace 


Para el prediseño del elemento de enlace se determina el cortante y 
momento último, en la sección del enlace, debido a la combinación de cargas 
descrita en 32, la misma se puede definir con las nomenclaturas de la ecuación 34 
para hacer referencia a Vy inksicmo Al Cortante en el elemento Link debido al sismo y 
Vi2-cm+cv hace referencia al cortante en el elemento Link debido a cargas 
gravitacionales. De manera similar se definen las variables para el momento flector 
en el elemento Link como M;ing y la ecuación 35 describe la combinación de 
estados de carga para determinar el momento flector. 


Viink = Viinksismo T V1.2-cm+cv (34) 
Vim =131+78 
Vine = 20.9T 
Miink = Miinksismo + Mi2-cm+cv (35) 


e 
ML inksismo = Viink 2 


(36) 
0.5 
Miink = (13.1 E =) + 3.5 


Miink = 6.7T:m 


Una vez obtenidos el cortante y momento flector actuantes en el elemento 
de enlace (Link), se procede al predimensionamiento de elementos por lo que se 
asume la sección IPE300 y se determina el momento y cortante plástico de la 
sección. 
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Datos de la sección 
IPE 300 
s 








Fy = 2530 kg/cm^2 


— 


Figura 15 Sección IPE300. 


Mp = Z- Fy = 602.1cm3 - 2531.052 21, 37 
A “> cm? 1000Kg 100cm (87) 


Mp = 15.24Tm 


Vp = 0.6 : Fy ' Aw 
p y (38) 
1cm i 1T 


Kg —— 
Vp = 0.6: 2531.05 — - (300 — 2: 10.7)mm- 7.1mm ( 1000 Kg 


10mm 


Vp = 30.04T 


En base a la longitud del enlace, se determina el factor para conocer si es 
un enlace puede fluir en corte, en flexión o una combinación de estas acciones. 


e 0.5 
Factor = ——— == 2— 
actor = Mp/Vp _ 15.24/30.04 (39) 
© __ 9.99 
Mp/Vp ` 


Dado que el Factor es 0.99 menor que 1.6 entonces la fluencia es por 
corte; el enlace se denomina como un elemento corto y considerando que la 
rotación máxima que puede alcanzar el enlace es 0.08 rad. La capacidad a corte 
se determina igual al corte plástico, y la capacidad a flexión es igual al momento 
plástico multiplicado por el factor @ = 0.90. En la tabla 7 se presenta un resumen 
de la demanda (D) y la capacidad (C) del enlace, también se presenta la relación 
D/C, la misma que es menor a 1. 


Tabla 7 Propiedades de la sección IPE300 























Mu = 6.7 Tm 
D-Mn = 13.7 Tm 
D/C Flexión = 0.49 
Vu = 20.9 T 
D-Vn = 27.0 T 
D/C Corte = 0.77 
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Una vez determinadas las características de demanda y capacidad se 
deben tomar en cuenta que el prediseño de la viga fuera del enlace, de los 
arriostramientos y las columnas se basa en la filosofía del diseño por capacidad, 
es decir, asumiendo que el elemento de enlace desarrollará la fluencia y disipación 
de energía mientras los otros elementos deben mantenerse esencialmente 
elásticos para este nivel de fuerzas (Figura 16). Las fuerzas debido a la fluencia 
del enlace sustituirán a la carga sísmica en las combinaciones de carga para el 
análisis. 

V 


t P 





Figura 16 Caso especial de análisis para el diseño por capacidad. 


Para cumplir con este principio la normativa recomienda diseñar la viga 
fuera del enlace para un cortante máximo probable en el enlace que se describe 
en la ecuación 40. La reducción del factor 0.88 se fundamenta en que la losa 
podría contribuir en la capacidad de la viga fuera del enlace. Para el diseño de las 
diagonales se determina el cortante máximo probable del enlace con la ecuación 
41. 


Vur = 0.88 - (1.25 - Ry- V, = 1.1: RY Vyk (40) 


Vur = 1.25 RY ` Vang (41) 


ña) 


Además, se asume que existirá una fuerza axial en el enlace Pyp, como se 
indica en la ecuación 42. 





Pup = 7 (42) 


Para el prediseño de las columnas se puede asumir de manera 
conservadora que todos los enlaces van a alcanzar el cortante máximo probable 
de la ecuación 40 (Figura 17), incluyendo el factor 0.88 debido a que es poco 
probable que todos los enlaces, de todos los niveles, ingresen al rango inelástico. 


Nuevas Funciones de CEINCI-LAB para análisis y diseño de pórticos ... 545 








LX 
Figura 17 Caso especial de análisis para el diseño por capacidad de las 
columnas. 


e Prediseño de la viga fuera del enlace 


El diagrama de cuerpo libre para el análisis de fuerzas en los elementos se 
muestra en la Figura 18. 


Vup21.10°R, Va 


lo L viga =| 


Pa Vyp"L/2H 


-| a/2 H 


- L/2 - 
Figura 18 Diagrama de cuerpo libre para el prediseño de la viga fuera del enlace. 


El momento flector (Memhyigq) Y la carga axial (Pemny;a) en la viga se 


pueden encontrar mediante las ecuaciones 48 y 49 respectivamente. El factor 
Fdyiga permite determinar la fracción del momento que debe ser resistido por la 
viga, las variables auxiliares se pueden determinar con las ecuaciones 43 a 47. 


E: Lex iga 
eu a rae (43) 
viga 
L-e 


ban (44) 
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Elgg 
Spiag = = (45) 
Diag 
Syiga 
Fdyiga = 2 — (46) 
dd Sviga t Spiag 
e 
Mentace = Vup 3 (47) 
M Emhyiga = Fadyiga * Mentace (48) 
V mp j L 
P Emhyiga = Pyp = 2-H (49) 


Donde Spiga Y Spiag SON las relaciones entre el módulo de elasticidad, la 
inercia del elemento respecto a la longitud de la viga y la diagonal 
respectivamente. Para los cálculos se asume la sección de viga IPE 500, cabe 
destacar que el elemento de enlace tiene dimensiones distintas, y para las 
diagonales se estima la sección cuadrada de acero con dimensiones 150mm x 
150mm x 10mm (Figura 19). 












































a s As TB 
Datos de la sección al a oar 
IPE 500 
s= 200 mm 
x X |, = 500 mm 
= 10.2 mm mo z 
t= 16 mm | | 
a k 
Fy = 2531.1 kg/cm^2 X a 
: phi = 0.9 B 
| Inercia = 1839 cm” 























(a) (b) 
Figura 19 a) IPE500; b) TB 150mm x 150mm x 10mm. 


En la tabla 8 se indica las solicitaciones y la capacidad a flexión, corte, 
axial y acciones combinadas flexión-axial, se corrobora que la sección IPE500 es 
adecuada para a viga fuera del enlace. 


Tabla 8 Demanda versus Capacidad de la viga fuera del enlace 






































Mu = Memnyiga + Mi2cm+cv = 44.1 Tm 
Ø:Mn = @-Mp = 48.0 Tm 
D/C Flexión = 0.92 
VP =Vyp = 43.0 T 
Ø:Vn =9:Vp = 65.2 T 
D/C Corte = 0.66 
Fcr = 24982 L 
Pu = Pēmhnyiga 7 41.8 T 
Ø: Pn = 251.2 T 
D/C Axial = 0.17 
D/C Flexión - Axial = 1.00 
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e  Prediseño de las diagonales 


En esta etapa para las diagonales se asume que la solicitación a flexión 
(Memhpiag) y axial (Pemhpiag) se determina con las ecuaciones 50 y 51, es 
importante notar que para las secciones impuestas de viga y diagonal, la inercia 
de la viga es muy superior a la de la diagonal y por ende la solicitación a flexión en 
la diagonal es pequeña. 


Spiag 
mate Sviga + Spiag 
Memhpiag = Fdpiag * Mentace (51) 
Pemhpiag = Pur * cos(0) + Vup ` sen(0) (52) 


En la tabla 9 se indica las solicitaciones y la capacidad a flexión, axial y 
acciones combinadas flexión-axial, se corrobora que la sección de la diagonal es 
adecuada. En esta etapa de prediseño la solicitación a flexión y axial debido a las 
cargas gravitacionales en la diagonal no se consideran. 


Tabla 9 Demanda versus Capacidad de la diagonal 





























Mu = Memnyiga + Mi2cm+cv = 1.2 Tm 
Ø:Mn = @-Mp = 6.7 Tm 
D/C Flexión = 0.18 
T 
Fcr = 17550 — 
m2 
Pu = Pemhyiga = 68.0 T 
Ø: Pn = 88.5 T 
D/C Axial = 0.77 
D/C Flexión - Axial = 0.93 














e Prediseño de las columnas 


Para prediseñar la columna más cargada ante cargas gravitaciones en 
planta baja, en este caso la columna del eje B, se calcula la carga axial en función 
del área colaborante de la columna de la Figura 10 y de las cargas gravitatorias 
antes definidas. La expresión 53 permite realizar la sumatoria de la carga de cada 
piso. 


num pisos 





| TOR . +L, 
Pu= > (1.2 -CMi+1.0- CVi) : (eripias _ Viga a viga derecta) (53) 


i=1 


El valor resultante es Pu = Pus »cm+cv = 154.66 (T). Para el estado de 
carga de sismo se considera la carga axial Pemn,,, CON la siguiente ecuación: 


num pisos 


P Emhco — 2 Vupi (54) 


i=1 
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Donde Vyp es el cortante máximo probable de los enlaces en cada piso, se 
determina su valor con la ecuación 52. 


Los valores de la demanda y capacidad de la columna se resumen en la 
tabla 10 (se asume una sección HEB 400 para la columna), para la capacidad 
axial se asume la relación de esbeltez de 60, se corrobora que la sección impuesta 
puede ser adecuada en la etapa de prediseño. 


Tabla 10 Demanda versus Capacidad de la columna 




















Fcr = 21057 L 
Pu = Puy zcm+cv + Pemhyiga = 282.3 T 
Ø-Pn = 403.3 T 
D/C Axial = 0.70 











A continuación, se verifica la relación entre momentos nominales 


probables de columna Y Mpc versus vigas } Mpv que convergen a una misma 
junta. 


> ce = 


Placas de Continuidad NM Articulación plástica 
p Mo. Articulación plástica 


M \ M 
J ~ 


V 


P, — 
Figura 20 Nudo de una columna interior de un pórtico. NEC-SE-AC 2015. 


En el articulo de Cagua et al. 2021 se demuestra relaciones aproximadas 
para cumplir con esta relación en base al módulo plástico de las secciones Z, para 
columnas este se entenderá como Zcol y para vigas como Zviga, y el esfuerzo de 
fluencia del material Fy. El nudo viga-columna de nivel 1, en el eje B, las 
secciones de los elementos se indican en la Figura 21 y se analizan de la 
siguiente manera. 


(55) 





i 
> e > 1.0 
Y Mpv 
Y 0.7: Zcol : Fy z 
y 1.15: Fy : Zviga — 





Y Zcol 
ST 2 1.64 
Y Zviga 
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HEB400 


o 1 
PEAS | E 500 | 


| HEB400 


Figura 21 Conexión viga-columna. 


La relación de la capacidad de las columnas versus la relación de 
capacidad de las vigas, en función de los módulos plásticos de sección de estos 
elementos, se indica en la tabla 11. En esta tabla se observa que esta relación es 
1.87 mayor a 1.64, por lo tanto, se asume que las secciones impuestas son 


adecuadas para proceder con el análisis. 


Tabla 11 Relación aproximada viga-columna nudo de análisis 























Zcol4gB4o0 = 21057 cm? 
Zvigđpesso = 1238.3 cm? 
ZvigđipEsoo = 2107.3 cm? 
Z_cols / Z_vigas = 1.87 ok 








El análisis se desarrollará con secciones similares a las presentadas en 
Cagua et al. 2021 y únicamente se modifica algunos elementos de vigas y 
diagonales como se muestra en la Figura 22, el pórtico corresponde a un sistema 
dual compuesto por elementos resistentes a momento y diagonales excéntricas. 


>N 


>X 4 2 3 


Figura 22 Secciones del pórtico con base en un prediseño. 


550 Brian Cagua, Roberto Aguiar, Julia Pilatasig y Alejandro Bonilla 





Las siguientes propiedades de material corresponden a las del acero 
ASTM A36, Fy = 36 ksi, Fu = 58 ksi y E = 29000 Ksi. Las dimensiones de las 
columnas y vigas se definen en la Tabla 12, donde la identificación del color es 
consistente con la Figura 23. 


Tabla 12 Descripción de las secciones de los elementos estructurales. 





Color Elemento Dimensiones (mm) 



































dxb x tf x tw 
Columnas 
' HEB400 400x300x24x13.5 
| Columnas 
tf HEB360 360x300x22.5x12.5 
Vigas 
IPE500 500x200x16.0x10.2 
Vigas 
IPE450 450x190x14.6x9.4 
` Vigas 
IPE400 400x180x13.5x8.6 
Vigas 
IPE300 300x150x10.7x7.1 
i Color Elemento Dimensiones (mm) 
T axe 
Diagonales 
HSS 150x10 
150X150X10 





5. MODELACION CON CEINCI-LAB 


Como se ejemplifica en Cagua et al, 2021 (1,2), para corroborar el 
prediseño se requiere un análisis numérico global de la estructura. Es necesario 
aplicar una metodología de análisis sísmico para determinar desplazamientos, 
derivas, fuerzas y cortante en cada piso; posteriormente verificar que las derivas 
sean menores a los límites permitidos por la NEC-15, para este caso debe ser 
menor a 2%; con las fuerzas en los pisos y con las cargas muertas y vivas se debe 
realizar un análisis estático para cada estado de carga, recordando que el efecto 
del sismo es bidireccional y con ello se obtienen las combinaciones de carga para 
calcular la envolvente de las solicitaciones (Demanda “D”). Posteriormente se 
debe determinar la capacidad axial, a corte, y flexión de los elementos, es 
importante conocer que las acciones a flexión y axial están interrelacionadas. Con 
estos valores se calcula la relación Demanda versus Capacidad “D/C” y se verifica 
que sea menor a 1. Además, en pórticos con arriostramientos excéntricos se debe 
realizar el diseño por capacidad, es decir reemplazar a los elementos de enlace 
con las fuerzas esperadas, que viene a ser el cortante máximo probable que 
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pueden llegar a desarrollan en función de la longitud del elemento. Estos efectos 
reemplazan al sismo, que en conjunto con las cargas muerta y viva en base a las 
combinaciones de carga, se debe verificar que la D/C sea menor a 1 para 
asegurar un comportamiento de enlace débil, viga y riostras fuertes y columna 
mucho más fuerte. Finalmente se verifica la conexión viga/columna. 


Los códigos del programa CEINCI-LAB para incluir diagonales con 
excentricidad en la viga central en el análisis han sido explicado en el trabajo de 
titulación de Romero (2020). En este artículo se presentará la metodología de 
análisis y se describirá con detalle las nuevas funciones de CEINCI-LAB para el 
diseño de pórticos con arriostramientos excéntricos. 


En color rojo de la Figura 23d se muestra la parte adicional del código para 
incorporar diagonales en V invertida en los marcos 2, 5, 8 y 11 (variable 
mar_Enlace), estos marcos son los presentados en la Figura 23c. Las rutinas 
dibujo, dibujogdl, dibujoNudElem, dibujovanos permiten graficar el esquema 
estructural; los grados de libertad en cada nudo y la numeración de nudos y 
elementos, además de la numeración de los vanos del pórtico en los que se puede 
colocar diagonales, esto se aprecia en las Figuras 23a, 23b y 23c. 


4 ‘ ‘ 4: 4 ‘ 


e + * . o + + ‘ 


(a) 


Esquema Estructural 
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10 1 12 
7 8 9 
4 5 6 
1 2 3 











(c) 





Código del programa principal, parte 1: 











rn GEOMETRÍA DEL PÓRTICO .........cccsssssssssseeseeeees 

sv =[5.00;7.00;5.00]; %Ingresar un vector con la longitud de vanos en (metros) 
Long_Piso=5*2+7.0; % Longitud (m) 

Ancho_Col=6; % Ancho colaborante 

sp =[3.6;2.88;2.88;2.88]; %Ingresar un vector con la altura de pisos en (metros) 
np = length(sp); % Numero de pisos 

Long_Enlace=0.5; % Longitud del elemento de enlace en metros 
mar_Enlace=[2;5;8;11];%Numero del marco con enlaces 
num_Enlaces=size(mar_Enlace,1); 

CM =[0.5;0.5;0.5;0.35]*Ancho_Col; % Carga muerta T/m en cada piso 
CV=[0.25;0.25;0.25;0.1]*Ancho_Col; % Carga viva T/m en cada piso 
Porc_viva=0.25; % Fracción de carga viva para el peso reactivo sísmico 
(W=D+Porc*L) 

%% Rutinas para generar la geometría del Pórtico 
[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_3nudos_viga(sv,sp,mar_Enlace); 
[X,Y]=glinea_portico_3nudos_viga(nv,np,sv,sp,nod,nr,Long_Enlace,mar_Enlace); 
[NI,NJ]=gn_portico_3nudos_viga(nr,nv,np,nudt,nudcol,nudvg,mar_Enlace); 

%% Incorporación de las diagonales 

num_mar=length(mar_Enlace); 

[NI2,NJ2]=geometria_3nudos_diagonales(nr, nv, np, nudt, mar_Enlace); 
[NI,NJ]J=gn_portico3(NI,NJ, Nl2,NJ2); 
[NI,NJ,Vigas_Interme]=rutina_paso_numeracion3(NI,NJ,nv,np,nudcol,nudvg, 
mar_Enlace); 

[CG,ngl]=cg_sismo2(nod,nr,Y); %Calcula los grados de libertad análisis sísmico 
for i=1:np 

qp(i,1)=-(CM(i,1)+(Porc_viva*CV(i,1))); % Carga distribuida en vigas en cada piso 
T/m 


end 
Voloradas laica DIBUJO DE LA ESTRUCTURA 2D..ooooccccnnncccccnnnnnnnnccos 
%dibujo(X,Y,NI,NJ) %Dibujo de estructura 


dibujogdl(X,Y,NI,NJ,CG)  %Dibujo de grados de libertad para análisis sísmico 
dibujoNudElem(X,Y,NI,NJ,CG) %Dibujo de estructura con numeración de nudos y 
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elementos 
dibujovanos(sv,sp); 





(a) 


Figura 23 Datos geométricos y restricciones del pórtico de análisis; a) Grados de 
libertad para análisis sísmico; b) Numeración de nudos y elementos; c) 
Numeración de vanos; d) Código de CEINCI-LAB con datos de 
geometría y cargas. 


Los tipos de secciones para columnas y vigas se asignan en la matriz 
Sec_VG_COL que contiene en cada fila un tipo de sección, además, en la primera 
columna de cada fila se indica el tipo de sección, seguido de las dimensiones del 
elemento; las dimensiones son peralte, d; espesor del alma, tw; ancho del patín 
superior, bs e inferior bi; espesor de los patines superior tfs e inferior tfi. 


Para designar cada sección a un elemento se debe tener en cuenta la 
numeración de estos. En la Figura 23b se muestra la numeración del pórtico. La 
matriz Elem_Sec_VG_COL permite asignar el tipo de sección y el elemento al que 
se le desea dar esta sección, colocando, en la primera columna el tipo de sección 
y en la columna siguiente el número de elemento al cual se asignara la misma. El 
código de CEINCI-LAB se expone en la Figura 24. 


Las secciones de diagonales se consideran en base a las propiedades de estas. 
Para calcular las propiedades de área, inercia y modulo plástico de la sección se 
utiliza la función Geom_acero que requiere de datos Sec_Col, que contiene: 
Material=1 para acero; Type=3 para sección tubular; d, tw, bs, bi, tfs, tfi las 
dimensiones del perfil; Lon es la longitud del elemento que se puede considerar 
como 3.0 aunque no es relevante en el cálculo; cpc es 1 si se trata de una 
conexión precalificada. Geom_acero calcula algunas propiedades de la sección 
que no se asignara a ninguna variable y por ello se tiene el símbolo -, mientras 
que, en los espacios correspondiente al Z, A, | se tienen variables. 


Las variables de interés se asignan a los elementos en base a 
Seccion_Diag, Z_Diag, Dimensiones _Diag, estas matrices se componen de 5 
columnas cada una; en la primera columna se coloca el número del elemento, al 
que se asignaran las propiedades en las columnas 2 y 3 de la matriz, (la 
numeración de los elementos se indica en la Figura 23b), la cuarta columna indica 
cuantos elementos adicionales tienen la misma propiedad y los siguientes 
elementos que secuencia de numeración tienen, por ejemplo se requiere asignar 
las propiedades tipo 1 al elemento 41 y adicionalmente a los elementos 42, 43 y 
44 (secuencia de 1 en 1). 





Código del programa principal, parte 2: 


iia SECCIONES: ¿cita 
% Propiedades del material 

Fy=25310.505; % Esfuerzo de fluencia, T/m2 
Es=20389019.16; % Módulo de elasticidad, T/m2 
Fu=40778.036; % Esfuerzo último, T/m2 
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Ry = 1.3; % Factor de Fluencia Probable (Ry) depende del acero para A36 
Ry=1.3 

%Uso Conexion = # % Conexion Precalificada RBS 1, Conexion 
Precalificada diferente que RBS 2, No Precalificada 3. 

Uso Conexion = 1; 

Materiales = [Fy, Fu, Es, Uso_Conexion]; 


% Sec_VG_COL = Elementos | =[tipo#, d,tw,bs,bi,tfs,tfi]... En milimetros!! 
% https://www.importaceros.com/ecuador-quito/viga-heb/ 


Sec_VG_COL=[ 1 400, 13.5, 300, 300, 24.0, 24.0; % Tipo 1 HEB 400 

2 360, 12.5, 300, 300, 22.5, 22.5; % Tipo 2 HEB 360 

3 450, 9.4, 190, 190, 14.6, 14.6; % Tipo 3 IPE 450 

4 500, 10.2, 200, 200, 16.0, 16.0; % Tipo 5 IPE 500 

5 400, 8.6, 180, 180, 13.5, 13.5; % Tipo 5 IPE 400 

6 300, 7.1, 150, 150, 10.7, 10.7; % Tipo 5 IPE 300 

7100, 1, 1, 1, 1, 13]; % Es para diagonales *No 
real la sección colocada, únicamente es para generar los elementos 


% Elem_Sec_VG_COL = ["tipo# respecto de Sec_VG_COL", Identificador de 
los elementos] 


Elem_Sec_VG_COL=[1 2, 6,10,14, 0, 0; 
1 3, 7,11,15, 0, 0; 
2 1, 5, 9,13, 0, 0; 
2 4, 8,12,16, 0, 0; 


317, 18, 21, 22, 23, 24; 
3 27, 28, 29, 30, 33, 34; 
419, 20, 25, 26, 31, 32; 
5 37, 38, 0, 0, 0, 0; 

5 35, 36, 39, 40, 0, 0; 

6 49, 50, 51, 52, 0, 0; 

7 41, 42, 43, 44, 45, 46; 
7 47, 48, 0, 0, 0, 0); 

Figura 24 Código de CEINCI-LAB con datos de materiales y secciones. 











En base a la función gelem_portico se construye una matriz de 2 
columnas, en las que se ordena las propiedades de interés de los elementos. 


En la Figura 25b se describe el código para ingresar las secciones de las 
riostras y en la variable Tconex se define que, si la riostra trabaja axialmente, este 
valor es cero; caso contrario estará sometida a efectos de flexión y axial. En 
pórtico con arriostramientos excéntricos se recomienda conexiones de las riostras 
a la viga que permitan transmitir parte del momento flector a estos elementos, 
como se indica en la Figura 25a. 
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Código del programa principal, parte 3: 








AAA SECCIÓN DIAGONALES eococncccncinnccnnonioncocnnooos 
% Diagonales tipo 1 TC 150x150x10 mm 
Material=1;Type=3;d=150/1000;tw=10/1000;bs=d;bi=bs;tfs=tw;tfi=tw; 
%Geometria de la Diagonal 


cpc=1; % Conexión precalificada 

Lon=3.0; % Longitud se requiere para calcular curvaturas 
Sec_Col=[Material,Fy ,Fu,Es, Type,d,tw,bs, bi, tfs, tfi]; 
[~.~,~,~ZA,L~,~,~,~,~,~,~]=Geom_acero(Sec_Col,Lon,cpc); 
A_Diagi=A; % Area de la seccion transversal 

|_Diagi =l; % Inercia de la seccion transversal 


h_Diag1=d-2*tw; % Altura de la diagonal para relación de aspecto 
tw_Diag1=tw; % Espesor del alma para relación de aspecto 
z_Diag1=Z; % Módulo plástico de la diagonal 


% Diagonales tipo 2 TC 150x150x10 mm 
Material=1;Type=3;d=150/1000;tw=10/1000;bs=d;bi=bs;tfs=tw;tfi=tw; 
%Geometria de la Diagonal 


cpc=1; % Conexión precalificada 

Lon=3.0; % Longitud se requiere para calcular curvaturas 
Sec_Col=[Material,Fy ,Fu,Es, Type,d,tw,bs, bi, tfs, tfi]; 

[~,~,~,~, ZAL~,~,~,~,~,~,~]=Geom_acero(Sec_Col,Lon,cpc); 
A_Diag2=A; % Area de la seccion transversal 

|_Diag2 = |; % Inercia de la seccion transversal 


h_Diag2=d-2*tw; % Altura de la diagonal para relación de aspecto 
tw_Diag2=tw; % Espesor del alma para relación de aspecto 
z_Diag2=Z; % Módulo plástico de la diagonal 


Seccion_Diag=[41 A_Diag1 |_Diag1 3 1; 
45 A_Diag2 |_Diag2 3 1]; % SECCION=Í[i, base, altura, nig,ii] 
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Z_Diag=[41 z_Diag1 03 1; 
45 z_Diag2 0 3 1]; % SECCION=T[i, z_Diag1, 0, nig,ii] 


Dimensiones_Diag=[41 h_Diag1 tw_Diag1 3 1; 
45 h_Diag2 tw_Diag2 3 1]; % Dimensiones_Diag=[i, h o Diametro, 
espesor, nig, li] 


[ELEM_DIAG]=gelem_portico(Seccion_Diag); % 
ELEM_DIAG=[Area_Diagonales,lnercia Diagonales] 

[ELEM_Z Diag]=gelem_portico(Z_Diag);  % Matriz con modulos plastico 
de diagonales 


Tconex=1; % Tipo de conexión de las diagonales (0 = articulada // 1=a 
momento) 


[num_Mar_Diag,variable]=size(mar_Enlace); % Número de marcos con 
diagonales 

for u=1:num_Mar_Diag 

num_viga(u,1)=nudcol+2*mar_Enlace(u,1)-1; 
num_viga(u,2)=nudcol+2*mar_Enlace(u, 1); 

end 











(b) 
Figura 25 Análisis del enlace en CEINCI-LAB; a) Diagrama de cuerpo libre en el 
enlace; b) Código de CEINCI-LAB con datos de materiales y 
secciones. 


Para el análisis sísmico se consideran los parámetros presentados en la 
Figura 26, en este caso la función 
Analisis_Sismico_Porticos Diag_Excentricas permite calcular las matrices de 
masa, rigidez, periodos de cada modo de vibración (modelo de 1 grado de libertad 
por piso), periodo con la ecuación de la NEC-15, cortante, fuerzas laterales, 
derivas y desplazamientos inelásticos de cada piso. La función se denomina 
Diag_Excentricas debido a que tiene diagonales con una excentricidad. 
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Código del programa principal, parte 5: 








Soo ica CALCULOS SISMICOS ...... eects 
% Espectro de diseño NEC 15 
is=4; % Suelo D 
iz=5; % Zona sísmica 5 = 0.4g 
ip=2; % Region Sierra 
Importancia = 1.0; % Factor de importancia 
R = 6; % Factor de reduccion de la fuerza sismica 
Cd = 5; % ASCE 7 - 16 SCBF Tabla 12.2-1 
omega=2; % Factor de Sobre resistencia SCBF( ASCE 7-16 Tabla 12.2-1 ) 
Dibujar = 1; % 1 para dibujar el espectro // 2 no dibujar 
fip=1 ;fie=1.0; % Factores de planta y elevacion 
zeda=0.05; % Amortiguamiento 
Tipo_Estr =2; %1 estructuras de acero sin arriostramiento; 2 con arriostramiento 
% 3 estructuras de hormigón sin muros ni diagonales; 4 con 
% muros y diadonales rigidizadoras y estructuras basadas en 
% estos 
% Tipo_Estr =1 Acero sin arriostramientos // Tipo_Estr =2 Acero con 
arriostramiento 
% Tipo_Estr =3 Hormigon armado sin arriostramientos // Tipo_Estr =4 Hormigon 
armado con arriostramiento 
Espectro=[is iz ip Importancia R fip fie zeda Tipo_Estr Dibujar]; % Datos para el 
espectro NEC 15 


aiii ANALISIS SISMICO DE LA ESTRUCTURA 

[M_Masa, M_Rigidez_Sis, T, T_Nec, Cort_Sis, Fuerzas_Sis, Derivas, 
Desp_Inelast]=Analisis_Sismico_Porticos_Diag_Excentricas(sv,sp,mar_Enlace,Lo 
ng_Enlace,qp,Sec_VG_COL,Elem_Sec_VG_COL,ELEM_DIAG,Materiales,Espectr 
o, Tconex); 


Vo orando PERIODOS saiia idad 
disp('Los periodos de vibración en segundos con CEINCI-LAB son:’) 
T % seg 


disp('El periodo de vibración fundamental en segundos NEC 15 es:') 
T_Nec % seg 





Figura 26 Código de CEINCI-LAB para análisis sísmico. 


La matriz de rigidez se muestra a continuación, se puede observar que la 
dimensión de esta matriz es de 4x4 debido a que se considera un modelo 
simplificado de 1 grado de libertad por piso, entonces la matriz de rigidez de la 
estructura con los grados de libertad indicados en la Figura 23a se condensa en 
los grados de libertad principales (los horizontales): 
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63825 -41460 4902.3 533.76 
Ve -41460 74878 -41583 4588.8 
4902.3 -41583 69559 -32295 

| 533.76 4588.8 -32295 27215 | 





La matriz de masa también tiene dimensión 4x4 debido a que el modelo 
concentra la masa del piso en un punto, es decir una masa en cada nivel. 





5.85 0 0 0 
iie 0 5.85 0 0 
0 0 5.85 0 

| oO 0 0 3.90 | 





El periodo fundamental de vibración calculado con la matriz de masa y 
rigidez es de 0.26 segundos, que es inferior al límite de 1.3 veces el periodo de la 
NEC 15. En la Figura 27 se puede visualizar el espectro elástico de diseño de 
color azul y de color verde el espectro reducido por el factor R, además la línea 
vertical de color rojo representa el periodo fundamental de vibración y las otras 
líneas de color gris son los otros periodos de vibración. 


La función de CEINCI-LAB permite calcular los desplazamientos y derivas 
inelásticas de cada piso, para la deriva se emplea la ecuación 56, considerando la 
variable Factor = 1, esto de manera conservadora para estructuras metálicas, 
también el usuario puede modificar la rutina si desea considera Factor = 0.75 
como lo indican algunas normativas. 


Am = Factor: Ae :R (56) 


_ Despl. Rel.) 


a (57) 


Ae 


Donde Am es la deriva inelastica para cada piso, este valor debe ser 
menor a 0.02 para pórticos de acero NEC-15 (2015); R es el factor de modificación 
de la respuesta sísmica; Ae es la deriva elástica para cada piso; Despl. Rel.) es 
el desplazamiento de cada piso obtenido con la aplicación de fuerzas laterales; h 
es la altura de cada piso. 


En la Figura 28 se muestra los resultados de desplazamientos, derivas y 
cortante en cada piso del pórtico, determinados con el Método Estático 
Equivalente en la línea de color rojo y la de color azul representa los resultados del 
Método Modal Espectral en base a la combinación cuadrática completa. Se puede 
apreciar que las respuestas del método Modal Espectral son inferiores a las 
respuestas del método Estático Equivalente. El usuario debe escoger con que 
fuerzas laterales continua el análisis, no obstante, se debe cuidar que la respuesta 
del análisis Modal Espectral sea superior al 80% del método Estático Equivalente 
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si el usuario decide emplearlo en los siguientes análisis. En este caso el cortante 
basal con el método Estático Equivalente es 41.7 T y con el método Modal 
Espectral es 37.64 T; esto representa 90 % entre un método y otro. 


DESPLAZAMIENTO DERIVA POR PISO CORTE POR PGO 


Figura 28 Respuestas del análisis sísmico. 


La deriva inelástica máxima que se presenta en el segundo piso es 0.3 % 
y 0.27 % con cada método, en ambos casos es menor al 2% que se tiene como 
límite en la NEC 15 para pórticos de acero. 


Posteriormente se debe calcular el índice de estabilidad de la estructura 
(0), esto se relaciona con el efecto P-Delta, es decir con efectos adicionales que 
pueden causar un incremento en las fuerzas internas, momentos y derivas, esto se 
debe considerar para la evaluación de la estabilidad estructural ante posible volteo 
y se determina con la ecuación 58. 





Donde Pi es la carga axial debido a los estados de cargas muerta y viva, 
es posible considerar la reducción de cargas vivas como se establece en las 
normativas, sin embargo, para ser conservadores en CEINCI-LAB no se considera 
reducción. A; es el desplazamiento relativo elástico de cada piso, si se considera 
una conexión RBS debido a la reducción de los patines en las vigas se puede 
aumentar un 10% estos desplazamientos como lo indica la NEC-SE-AC 2015; V; 
es el cortante en cada nivel; h; es la altura de cada nivel. En la Figura 29 se 
presenta el código de CEINCI-LAB, la función Indice _Estabilidad_Sismo_V2 
utiliza para el cálculo de la estabilidad las fuerzas y desplazamientos obtenidas 
con el método Estático Equivalente, esto se define en la variable Metodo que es 
igual a 2 en este caso, si se requiere calcular con el método modal espectral esta 
variable debería definirse como 1. 
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Código del programa principal, parte 6: 





o aia ÍNDICE DE ESTABILIDAD ....oooococccncncncncncninnaciccnnncnss 
Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Lineal Estático Nec 15 


[Teta, Teta2, Teta_Admisible]=Indice_Estabilidad_Sismo_V2(CM, CV, 
Long_Piso, sp, R, Metodo, Cort_Sis, Desp_Inelast, Uso_Conexion); %NEC 15 











Figura 29 Código para análisis de estabilidad con CEINCI-LAB. 


Cuando el índice de estabilidad es menor a 0.10 entonces se clasifican los 
pisos como no susceptibles al volteo; si este índice es superior a 0.3 entonces la 
estructura es potencialmente inestable y debe rigidizarse; cuando se encuentra 
entre 0.1 y 0.3 se debe calcular un factor fp_, y multiplicar todas las fuerzas 
laterales por este factor. 


1 


fe-a = 14 (59) 


Siendo @ el mayor valor de los 6;. A continuación, se muestran los 
resultados del índice de estabilidad de cada piso de la estructura y se comprueba 
que es menor a 0.10 en todos los niveles, además, debido a la gran rigidez lateral 
de la estructura en comparación con pórticos a momentos entonces estos valores 
son pequeños. o o 

0.0049 
0.0040 
0.0026 


0.0013 








En esta etapa se puede indicar que las secciones de columnas, vigas y 
diagonales permiten controlar los desplazamientos laterales de manera adecuada, 
entonces se procede a calcular la relación de los módulos plásticos de columnas 
respecto a vigas que llegan a un mismo nudo, asegurando que sea mayor a 1.64 
como se demostró anteriormente. La rutina que calcula esta relación se denomina 
relacion_sum_Z columna_viga y dibujo relacion_sum_Z columna_viga 
permite presentar gráficamente esta relación (Figura 30), es importante mencionar 
que los códigos sísmicos indican que en el último piso no se requiere cumplir con 
esta relación de la capacidad de las columnas superior a la capacidad de las 
vigas. 








Código del programa principal, parte 7: 








%% Índice de Columna-Fuerte viga Debil 

rta cid 

[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y coseno 
num_Diag_Total =size(mar_Enlace,1)*2; 
nuVG_COL_DIAG=nudt+num_Diag_Total; 

% Inicio y fin de cada tipo de elemento 
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Inic_Fin_COL_VG=[1,nudt-num_Enlaces]; 
Inic_Fin_DIAG=[Inic_Fin_COL_VG(1,2)+1,Inic_Fin_COL_VG(1,2)+num_Diag_Tota 
I]; 
Inic_Fin_Enlace=[Inic_Fin_DIAG(1,2)+1,nuVG_COL_DIAG]; 
[Seccion,Elem_Tipo_Prop2,Z_ELEM,th_ELEM]=Diferentes_Vg_ COL(nuVG_COL_ 
DIAG,Sec_VG_COL,Elem_Sec_VG_COL,Uso_Conexion,Fy,Es,L); 
%Seccion=[ELEMENTO, Area, Inercia, NUMEROS DE ELEMENTOS A CREAR, 
INCREMETO EN NUMERO DE ELEMENTO] 

[ELEM]=gelem_portico(Seccion); 





Capacidad_sum_Z_columna_viga=relacion_sum_Z_columna_viga(coseno,nr,nudc 
ol,X,Y,NI,NJ,Z_ELEM,Elem_Tipo_Prop2,Sec_VG_COL); 
dibujo_relacion_sum_Z_columna_viga(X,Y,NI,NJ,Capacidad_sum_Z_columna_vig 
a) 





(a) 





Colg = 2.1 Coig = 1.3 ColVg = 1.3 ColWVg = 2.1 
ColVg = 3,2 Coig = 1.7 ColVg = 1.7 ColVg = 3.2 
Cog = 3,2 Caig = 1,7 Colg = 1,7 Colg = 3.2 
Coig =3.2 Cog = 1.7 Cog = 1,7 ColVg = 3.2 


(b) 
Figura 30 a) Código de CEINCI-LAB; b) Relaciones de módulos plásticos de 
columnas para vigas de forma gráfica. 


En esta etapa del análisis se comprueba que las secciones de columnas y 
vigas asignadas al pórtico dual permiten cumplir con los requisitos “viga débil y 
columna fuerte”, estabilidad global y derivas, por lo tanto, se puede continuar con 
el análisis estático, es decir aplicar los estados de carga gravitacionales y laterales 
(sismo) para determinar las fuerzas internas en los elementos. 


Para el análisis estático se cambian los grados de libertad y se considera 
un modelo con 3 grados de libertad por nudo (Desplazamiento en horizontal, 
vertical y Giro), esto debido a que es un pórtico plano. En la Figura 31a se muestra 
el código para calcular los grados de libertad con la rutina cg, es importante notar 
que se requiere definir las restricciones en la base de cada columna, para este 
caso se considera un empotramiento. En la Figura 31b se muestra la numeración 
de grados de libertad con la función dibujogdl. Con este modelo se calcula la 
matriz de rigidez KTT que es la suma de la matriz de rigidez del pórtico K_ACP 
(columnas y vigas) más la matriz rigidez de las diagonales K_ACD y el aporte de 
la rigidez de los elementos de enlace K_ACE. 
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Código del programa principal, parte 8: 








% Restricciones de apoyo 
aux=length(sv); 
for i=1:aux+1 
RES(i,:)=[i 1 1 1]; %Nudos restringidos del pórtico con diagonales 
end 
[CG,ngl]=cg(nod,nr,RES); %Calcula los grados de libertad 3 por nudo 
dibujogdl(X,Y,NI,NJ,CG) %Dibujo de grados de libertad para análisis estático 


[VC]=vc(NI,NJ,CG); %Vector de colocación 


%% Matriz de Rigidez 
% Contribución de Elementos de Vigas y Columnas 
ELEM_ACP=ELEM(Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2),:); 
L_ACP=L(:,Inic_Fin_COL_VG(1,1):Inic_Fin COL_VG(1,2)); 
seno_ACP=seno(:,Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2)); 
coseno_ACP=coseno(:,Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2)); 
% ... En la union de la viga - columna se libera momentos para hacer una 
% conexion a corte que permite grandes rotaciones 
VC_ACP=VC(Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2),:); 
iris 
[K_ACP]=krigidez_acero(ngl,ELEM_ACP,L_ACP,seno_ACP,coseno_ACP,VC_AC 
P,Es); %Matriz de rigidez de columnas y vigas 
% Contribución de Elementos Diagonales de Acero 
L_ACD=L(:,Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2)); 
seno_ACD=seno(:,Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2)); 
coseno_ACD=coseno(:,Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2)); 

if Tconex==0 

ELEM_ACD=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),1); 


TEE IMPORTANTE Solo Area 
VC_ACD=VC(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),[1,2,4,5]); 
ici IMPORTANTE IGNORO GRADOS DE LIBERTAD DE GIRO 


[K_ACD]=krigidez(ngl, ELEM_ACD,L_ACD,seno_ACD,coseno_ACD,VC_ACD,Es); 


VC_ACD2=[VC(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),[1,2]),0*VC(Inic_Fin_DIAG 
(1,1):Inic_Fin_DIAG(1,2),[3]),VC(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),[4,5]),0*V 
C(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),[6])]; 
else 
ELEM_ACD=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 
A e IMPORTANTE Solo Area 
VC_ACD=VC(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 
lola IMPORTANTE TODOS LOS GRADOS DE LIBERTAD 


[K_ACD]=krigidez_acero(ngl,ELEM_ACD,L_ACD,seno_ACD,coseno_ACD,VC_AC 
D,Es); 

VC_ACD2=VC_ACD; 

end 
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% Contribución de Elementos de Enlace 
[ELEME]=ELEM(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:)¡nmiembro=Inic_Fin_ 
Enlace(1,1); 
LE=L(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2));senE=seno(Inic_Fin_Enlace(1,1): 
Inic_Fin_Enlace(1,2));cosE=coseno(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1 ,2));VC 
E=VC(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:); 
[K_ACE]=krigidez_arm(ngl,ELEME,LE,senE,cosE,VCE,Es,CG,NI,NJ,nr,nmiembro) 


, 


KTT=K_ACP+K_ACD+K_ACE; %<------- MATRIZ "K" DE TODA LA 
ESTRUCTURA (T/m) 





(a) 


Grados de Libertad 
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(b) 
Figura 31 a) Código de CEINCI-LAB; b) Grados de libertad para el análisis 
estático. 


A continuación, se aplica las cargas sísmicas en la estructura, como se indica en 
la Figura 32, recordando que un caso de carga es cuando se analiza con fuerzas 
laterales de izquierda a derecha y otro caso es aplicando las fuerzas en sentido 
opuesto, esto es importante en estructuras que no son simétricas. 


F4 F4 


— — 
FS 4" 
El dl «2 

a Fl 


Figura 32 Cargas laterales aplicadas a las estructuras. 
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En la Figura 33 se muestra el código de CEINCI-LAB que permite realizar 
el análisis estático con fuerzas laterales, el procedimiento consiste en definir las 
fuerzas en los nudos, esto se lo hace en el vector F. Se calcula el vector de cargas 
generalizadas, es decir las cargas en los nudos de las estructuras que son 
correspondientes a los grados de libertad indicados en la Figura 31b. Se resuelve 
la operación matricial planteada en la ecuación 60, donde Q y K son el vector de 
cargas generalizadas y la matriz de rigidez de la estructura (Aguiar 2020). 
Entonces se calculan los deslazamientos y giros del vector q (se debe considerar 
los grados de libertad de la Figura 31b). 


Q=K-q (60) 


Cuando la estructura tiene riostras y elemento de enlace, se debe dividir 
en el análisis para el pórtico con columnas y vigas, el análisis para las diagonales 
(las rutinas _Articuladas hacen referencias a la liberación de giros en las vigas) y 
posteriormente los elementos de enlace. 








Código del programa principal, parte 9: 








nds CARGAS SÍSMICAS .....c:cccsssscsessessesesssseeees 
Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Lineal Estático Nec 15 
njc=np; % Número de juntas Cargadas 
% 
Ypisos=nonzeros(unique(Y)); 
oij=1; 
uij=1; 
for i=1:nod 
if X()==0 88 any(Y(i)==Ypisos) 
Nudos_Izquierda(oij)=1; 
oij=1+o0ij; 
end 
if X(i)==max(X) 88 any(Y(i)==Ypisos) 
Nudos_Derecha(uij)=i; 
uij=1+uij; 
end 
end 
%% Sentido de izquierda a derecha (Positivo) 
F = [[Nudos_Izquierda]' Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % F = 
[Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=0; % Número de miembros cargados 
Fm=0; %Elem carg, carga, código, elem a gener, incr numero elemento 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ);% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Sismo=KTTIQ_Sismo; 
% Fuerzas y momentos finales en los elementos 
nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
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[Fuerzas Elementos _Sismo_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELEM_AC 
P,L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Sismo,Q2_Sismo(Inic_Fin_COL 
_VG(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 


nmiembro= 1; % Numero del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

ELEM_ACD2=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 

RES_ACD=0; 

% Fuerzas y momentos finales en las diagonales 


[Fuerzas_Elementos_Sismo_DG]=fuerzas_acero_arm1_V2(ngl, ELEM_ACD2,L 
_ACD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Sismo,Q2_Sismo(Inic_Fin_DIAG( 
1,1):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 


% Fuerzas en enlace 
nmiembro=1; 
RES_E=0; 


[FA_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cosE,VCE,Es,q_Sismo,Q 
2 Sismo(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ,nmiembro,RES 
ER 

% Fuerzas y momentos finales en toda la estructura 

Fuerzas Elementos_Sismo_P=[Fuerzas_Elementos_Sismo_COL_VG;Fuerzas_ 
Elementos_Sismo_DG;FA_EN]; 








Figura 33 Código de CEINCI-LAB. 


Con el vector de coordenadas generalizadas q y la matriz de rigidez de 
cada elemento k se calcula en coordenadas locales, las fuerzas en cada 
elemento; cada elemento tiene 6 grados de libertad, 3 en el nudo inicial y 3 en el 
nudo final. La matriz de celdas Fuerzas Elementos Sismo contiene en su 
posición {1,1} las fuerzas en los elementos debido a las fuerzas laterales de 
izquierda a derecha y en la posición (2,1) las fuerzas en los elementos debido a 
fuerzas laterales en sentido opuesto como se muestra en la Figura 34. 





Código del programa principal, parte 10: 











%% Sentido de derecha a izquierda (Negativo) 
F = [[Nudos_Derecha]' -Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % F = 
[Nudo cargado, FH, FV, Momento] 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ);% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Sismo=KTTIQ_Sismo; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
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[Fuerzas_Elementos_Sismo_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACP, 
L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Sismo,Q2_Sismo(Inic_Fin_COL_V 
G(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,amiembro,RES); 


nmiembro= 1; % Numero del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

ELEM_ACD2=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 

RES_ACD=0; 

% Fuerzas y momentos finales en las diagonales 


[Fuerzas_Elementos_Sismo_DG]=fuerzas_acero_arm1_V2(ngl, ELEM_ACD2,L_A 
CD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Sismo,Q2_Sismo(Inic_Fin_DIAG(1,1 
):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 

% Fuerzas en enlace 

nmiembro=1; 

RES_E=0; 


[FA_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cosE,VCE,Es,q_Sismo,Q2 
_Sismo(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ,nmiembro,RES _ 
E); 

% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_Sismo_N=[Fuerzas_Elementos_Sismo_COL_VG;Fuerzas_ El 
ementos_Sismo_DG;FA_EN]; 

% Fuerzas y momentos finales sismic@s positivo=P y negativo=N en la estructura 
con diagonales 

Fuerzas_Elementos_Sismo{1,1}=Fuerzas_Elementos_Sismo_P; 
Fuerzas_Elementos_Sismo{2,1}=Fuerzas_Elementos_Sismo_N; 








Figura 34 Código de CEINCI-LAB. 

Para el caso de cargas muertas y vivas se realizan análisis individuales de 
estados de carga como se indica en la Figura 35, se reitera que el problema se 
divide en encontrar las fuerzas en los elementos del pórtico principal, 
posteriormente calcular las fuerzas en los elementos diagonales y finalmente en 
los elementos de enlace e integrarlos al final del análisis en las matrices de 
Fuerzas_Elementos_CM y Fuerzas_Elementos_CV. 





Código del programa principal, parte 11: 











RO CARGAS GRAVITACIONALES „osii 
% Identificar en el piso donde se encuentra cada diagonal 
mar_pis=nr-1; % vanos en cada piso 
for j=1:np 
for i=1:num_Enlaces 
if j== 
if mar_Enlace(i)<=mar_pis* 
Piso_Enlaces(i, 1)=j; 
end 
else 
if mar_Enlace(i)>mar_pis*(j-1) && mar_Enlace(i)<=mar_pis*j 
Piso_Enlaces(i,1)=j; 
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end 
end 
end 
end 


njc=0; % Número de juntas Cargadas 
F = 0; % F = [Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=nudvg; % Número de miembros cargados 
ui=1; 
for i=1:np 
for j=1:nv*2 
Fm_CM(ui,:)=[nudcol+ui CM(i,1) 1 0 0]; %Elem carg, carga, código, elem a 
gener, incr numero elemento 
Fm_CV(ui,:)=[nudcol+ui CV(i,1) 1 O 0]; %Elem carg, carga, código, elem a 
gener, incr numero elemento 
ui=ui+1; 
end 
end 
for u=1:num_Enlaces 
Fm_CM(ui,:)=[Inic_Fin_Enlace(1,1)+u-1 CM(Piso_Enlaces(u,1),1) 1 0 0]; %Elem 
carg, carga, código, elem a gener, incr numero elemento 
Fm_CV(ui,:)=[Inic_Fin_Enlace(1,1)+u-1 CV(Piso_Enlaces(u,1),1) 1 0 0]; %Elem 
carg, carga, código, elem a gener, incr numero elemento 
ui=ui+1 ; 
end 


[Q_CM,Q2_CM]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm_CM,datos, 
nr,NI,NJ); 
[Q_CV,Q2_CV]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm_CV,datos,nr 
sNI,NJ); 

% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 

% Desplazamientos y Giros 

q_CM=KTT\Q_CM; 

q_CV=KTT\Q_CV; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Numero del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 


%% Fuerzas y momentos finales por carga muerta == CM 

% Fuerzas y momentos finales en columnas y vigas 
[Fuerzas_Elementos_CM_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACP,L_ 
ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_CM,Q2 _CM(Inic_Fin_COL_VG(1,1):l 
nic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 

% Fuerzas y momentos finales en diagonales 
[Fuerzas_Elementos_CM_DG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACD2,L_ACD, 
seno_ACD,coseno_ACD,VC_ACD2,Es,q_CM,Q2_CM(Inic_Fin_DIAG(1,1):Inic_Fin 
_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 

% Fuerzas y momentos finales en Enlaces 
[Fuerzas_Elementos_CM_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cos 
E,VCE,Es,q_CM,Q2_CM(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ 
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snmiembro,RES_E); 

% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_CM=[Fuerzas_Elementos_CM_COL_VG;Fuerzas_Elementos 
_CM_DG;Fuerzas_Elementos_CM_EN]; 


%% Fuerzas y momentos finales por carga viva == CV 

% Fuerzas y momentos finales en columnas y vigas 
[Fuerzas_Elementos_CV_COL_VGj]=fuerzas_acero_arm1_V2(ngl,ELEM_ACP,L_ 
ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_CV,Q2_CV(Inic_Fin_COL_VG(1,1):1 
nic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 

% Fuerzas y momentos finales en diagonales 
[Fuerzas_Elementos_CV_DG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACD2,L_ACD, 
seno_ACD,coseno_ACD,VC_ACD2,Es,q_CV,Q2_CV(Inic_Fin_DIAG(1,1):Inic_Fin 
_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 

% Fuerzas y momentos finales en Enlaces 
[Fuerzas_Elementos_CV_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cosE 
,VCE,Es,qg_CV,Q2_CV(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ,n 
miembro,RES_E); 

% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_CV=[Fuerzas_Elementos_CV_COL_VG;Fuerzas_Elementos 
_CV_DG;Fuerzas_Elementos_CV_EN]; 














Figura 35 Código de CEINCI-LAB. 


Basado en el principio de superposición los efectos de cargas en los 
elementos se suman directamente entre los diferentes estados de carga para tener 
las combinaciones de carga. En la NEC-SE-CG 2015 se tienen las siguientes 
combinaciones de carga al considerar el método LRFD: 


Combinación 1 1.4D 

Combinación 2 1.2D + 1.6L + 0.5max[Lr; S; R] 
Combinación 3 1.2D + 1.6max[L; S; R] + max[Lr; 0.5W] 
Combinación 4 1.2D + 1.0W + L + max[Lr; S; R] 
Combinación 5 1.2D +1.0E+L+0.25 
Combinación 6 0.9D + 1.0W 
Combinación 7 0.9D + 1.0E 


Donde D corresponde a la carga muerta (carga permanente + peso 
propio); E carga de sismo (componente horizontal y vertical); L sobrecarga viva; Lr 
es la sobrecarga viva en cubierta; S carga de granizo; R es la carga de lluvia; W 
carga de viento. 


En las combinaciones de carga donde se incluye los efectos del sismo es 
posible hacer una reducción del efecto de la carga viva en función de la 
probabilidad de ocupación, no obstante, no se considera esta reducción para ser 
conservadores en el análisis. También se debe considerar que las combinaciones 
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de carga con sismo implican, la componente vertical, pero en estos cálculos no se 
incluye este efecto. Con estas consideraciones se calcula las combinaciones de 
carga aplicables en este caso particular, con los efectos de carga muerta, viva y 
sismo (con fuerzas laterales en ambos sentidos); en la matriz Combinaciones se 
tiene las fuerzas en los elementos con las combinaciones de estados de carga y 
en la matriz Envolvente se determina los máximos efectos de estas 
combinaciones en cada elemento (máximos y mínimos). El código de CEINCI-LAB 
se presenta en la Figura 36, se indica que la función se denomina 
Combinacion_Cargas. 








Código del programa principal, parte 12: 








[Combinaciones, 
Envolvente]=Combinacion_Cargas(Fuerzas_Elementos_CM,Fuerzas_Elementos 
_CV,Fuerzas_Elementos_Sismo); 





Figura 36 Código de CEINCI-LAB. 


Los efectos del sismo con el factor de la sobre resistencia se consideran 
con el estado de carga indicad en la ecuación 61. 


Emn = 2° Qe (61) 


Donde Emn es la componente horizontal del sismo con el efecto de sobre 
resistencia; Q es el factor de sobre resistencia (en el ASCE 7-16 se exponen 
valores referenciales en función del sistema estructural); Qg es el efecto de sismo 
horizontal (cortante basal). 


Asumiendo que se desprecia el efecto de sismo vertical entonces las 
combinaciones de carga 5 y 7 cambian como se indica a continuación. 


Combinación 5b 1.2D + 1.0E mn + L + 0.2S 
Combinación 7b 0.9D + 1.0E mn 


Para considerar estos nuevos estados de carga se hace un análisis similar 
al caso anterior, con la particularidad de que se debe mayorar las fuerzas laterales 
debido al sismo en su componente horizontal, el código que genera esto se 
muestra en la Figura 37. 








Código del programa principal, parte 13: 








Volo socia CARGAS SÍSMICAS CON FACTOR DE 
SOBRERESISTENGIA sissiscetseasiacertoaidsecastanscare 

Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Lineal Estático Nec 15 
njc=np; % Numero de juntas Cargadas 
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%% Sentido de izquierda a derecha 
F = [[Nudos_lzquierda]' omega*Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % 
F = [Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=0; % Numero de miembros cargados 
Fm=0; %Elem carg, carga, código, elem a gener, incr numero elemento 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 
% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Sismo_omega=KTTIQ_Sismo; 
nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

% Fuerzas y momentos finales en los elementos 


[Fuerzas_Elementos_Sismo_omega_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELE 
M_ACP,L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Sismo_omega,Q2_Sismo( 
Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 

% Fuerzas y momentos finales en las diagonales 





[Fuerzas_Elementos_Sismo_omega_DG]=fuerzas_acero_arm1_V2(ngl,ELEM_AC 
D2,L_ACD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Sismo_omega,Q2_Sismo(Ini 
c_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 

% Fuerzas en enlace 





[Fuerzas_Elementos_Sismo_omega_EN]=fuerzas_acero_arm1_V2(ngl, ELEME,LE 
,senE,cosE,VCE,Es,q_Sismo_omega,Q2_Sismo(Inic_Fin_Enlace(1,1):Inic_Fin_Enl 
ace(1,2),:),CG,nr,NI,NJ,nmiembro,RES_E); 

% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_Sismo_P_omega=[Fuerzas_Elementos_Sismo_omega_COL 
_VG;Fuerzas_Elementos_Sismo_omega_DG;Fuerzas Elementos _Sismo_omega 
EN]; 














%% Sentido de derecha a izquierda 
F = [[Nudos_Derecha]' -omega*Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % 
F = [Nudo cargado, FH, FV, Momento] 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 
% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Sismo_omega=KTTIQ_Sismo; 
nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

% Fuerzas y momentos finales en los elementos 


[Fuerzas_Elementos_Sismo_N_omega_COL_VG]=fuerzas_acero_arm1_V2(ngl,E 

LEM_ACP,L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Sismo_omega,Q2_Sis 

mo(Inic_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 
% Fuerzas y momentos finales en las diagonales 





[Fuerzas_Elementos_Sismo_N_ omega _DG]=fuerzas_acero_arm1_V2(ngl,ELEM_ 
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ACD2,L_ACD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Sismo_omega,Q2_Sismo( 
Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ ACD); 
% Fuerzas en enlace 


[Fuerzas_Elementos_Sismo_N_omega_EN]=fuerzas_acero_arm1_V2(ngl,ELEME, 
LE,senE,cosE,VCE,Es,q_Sismo_omega,Q2_Sismo(Inic_Fin_Enlace(1,1):Inic_Fin_ 
Enlace(1,2),:),CG,nr,NI,NJ,nmiembro,RES_E); 

% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_Sismo_N_omega=[Fuerzas_Elementos_Sismo_N_omega_C 
OL_VG;Fuerzas_Elementos_Sismo_N_omega_DG;Fuerzas_Elementos_Sismo_N 
_omega_EN]; 











% Fuerzas y momentos finales sismic@s positivo=P y negativo=N en la estructura 
con diagonales 
Fuerzas_Elementos_Sismo_omega{1,1}=Fuerzas_Elementos_Sismo_P_omega; 
Fuerzas_Elementos_Sismo_omega{2,1}=Fuerzas_Elementos_Sismo_N_omega; 


VO taa COMBINACIONES DE CARGA CON 
SOBRERRESISTENCIA eer 

[Combinaciones_omega, 
Envolvente_omega]=Combinacion_Cargas(Fuerzas_Elementos_CM,Fuerzas_Ele 
mentos_CV,Fuerzas_Elementos_Sismo_omega); 





Figura 37 Código de CEINCI-LAB. 


En la sección 6.2 de la NEC-SE-AC 2015 se manifiesta que columnas, 
vigas y diagonales de pórticos especiales deben tener secciones clasificadas 
como sismicamente compactas (para pórticos duales), esto se puede interpretar 
como secciones de alta ductilidad como lo indica el AISC 341-16 en su capítulo D. 
Este requerimiento es importante debido a que asegura que las secciones no sean 
susceptibles de fallas anticipadas por inestabilidad local o pandeo local. 


En la tabla D1.1. del AISC 341-16 se muestran los límites para clasificar a 
las secciones según la relación de aspecto de patines y alma denominada 4; si 
este factor es menor a Ana la sección es altamente dúctil y si se encuentra entre 
Ana Y Ama se clasifica como moderadamente dúctil; si el valor de A es superior a 
este límite entonces la sección es no compacta. 
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Tabla 11 Limites de ancho versus espesor para elementos sometidos a 
compresión. Fragmento de la tabla D1.1 del AISC 341-16. 
Límites de 
A Referencia 
Ana Ama 
Elemento de Alta Elemento de Moderada 
ductilidad ductilidad 
pS 3 
TÍ 
E E 0.32 sl 0.40 = 
a| t Ry: Fy Ry: Fy 
Para Ca < 0.114 Para Ca < 0.114 
2.57 TAT (1 — 1.04C4) 3.96 Ry" Fy — 3.04C,) 
Para Ca > 0.114 Para C, > 0.114 
0.88 a 2.68 —C 1.29 sa 2.12 -C 
z h j Pe ~ a) . Ry Fy a a) 
< | tw ee: sag E 
> 1.57 Ry Fy > 1.57 D Fy 
Donde Donde 
C= Pu C= Pu 
Øe: Py “Oc * Py 
Py = Ry: Fy: Ag Py = Ry : Fy: Ag 
YN 
2 bo. 
> b 0.65 7 0.76 sl E 
2 t i Ry: Fy i Ry: Fy 
2 
Oo 
a 
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El Cálculo se realiza con la función Compacidad_Elementos_Diag 
indicada en la Figura 38a, como resultados se tiene la matriz compacidad que 
almacena la clasificación del patín, del alma y de la sección de cada elemento con 
códigos: 1 para HD, 2 para MD y O para ND, en este caso HD significa que la 
sección es altamente dúctil y se tuviera la nomenclatura MD es moderadamente 
dúctil, finalmente para secciones con poca ductilidad se establece la nomenclatura 
ND. La matriz Lamba contiene los límites para HD o MD de patín y alma de cada 
elemento. Finalmente, en la matriz denominada relacion_aspecto se almacena la 
relación de ancho para espesor de patín y alma de cada elemento. Además, se 
presenta gráficamente los resultados con la función dibujo_Compacidad como se 
muestra en la Figura 38b, los códigos HD, MD y ND también definen un color al 
elemento, como azul, verde y rojo respectivamente. Para este caso todas las 
secciones son altamente dúctiles. 








Código del programa principal, parte 14: 








YY ..... COMPACIDAD DE ELEMENTOS CONSIDERANDO SMF * AISC -341-16 
ELEM_COMP=[ELEM_ACP;ELEM_ACD2;ELEME]; 

[Compacidad, Lambda, 
relacion_aspecto]=Compacidad_Elementos_Diag_Excentricas(nudcol,Inic_Fin_CO 
L_VG,Inic_Fin_DIAG,Inic_Fin_Enlace,Es,Fy,Ry,ELEM_COMP,Elem_Tipo_Prop2,S 
ec_VG_COL,Envolvente_omega,Dimensiones_Diag); 
dibujo_Compacidad(X,Y,NI,NJ, Compacidad) 

% HD=Altamente ductil // MD = Moderadamente Ductil // ND = No Ductil 





(a) 
Esquema Estructural: Compacidad 
HO HD HO -HO +D HO HD 
HO HD AD HD HO D 
Ho HD HO7H0—HD—|—HD' #0- 
| 
HD HD HD “o Ho H 
| 
o HB HD- HBD -HD oD i HD mo 
HD HD AD HD HD HD 
o HD mo HO -HD PD t HB mo 
HO HD AD eo i HD 
(b) 
Figura 38 a) Código de CEINCI-LAB; b) Clasificación de las secciones en el 


pórtico. 


Para calcular la capacidad axial de los materiales se requiere del factor de 
longitud efectiva K, el mismo depende de las condiciones de borde del elemento. 
En el artículo de Cagua et al, 2021 (1 y 2) se ejemplifica el cálculo manual de este 
parámetro para los pórticos a momento y con arriostramientos concéntricos. 
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Para las vigas se considera un factor K=1 y para diagonales K=0.75. En la 
Figura 39 se muestra el código para calcular la capacidad de los elementos y para 
comparar la demanda versus la capacidad axial, a corte y flexión con las diferentes 
combinaciones de carga (el procedimiento se detalla en Cagua et al, 2021). Se 
destaca que únicamente cambian las rutinas Factor_K_Long_Efectiva_Diag y 
Capacidad_Elementos_Acero Diag_Excentricas, los cambios se dan para 
considerar diagonales en el análisis. 








Código del programa principal, parte 15: 








So Vo roriccaa CAPACIDAD DE ELEMENTOS DE ACERO ..........:::ececeeeeeeeeeeeeeeeteeeees 
Z_ELEM_COMP=[Z_ELEM]; 
Z_ELEM_COMP(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2))=ELEM_Z_Diag(Inic_Fin 
_DIAG(1,1):Inic_Fin_DIAG(1,2),1); 

for j=1:num_Diag_Total 

Seccion_Diag2(j, 1:5)=[Inic_Fin_DIAG(1,1)+j-1 ELEM_DIAG(Inic_Fin_DIAG(1,1)+j- 
1,1:2) j 0]; % Nota en esta matriz la 4ta columna numera las diagonales 

end 

Seccion2=[Seccion]; 
Seccion2(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),1:5)=Seccion_Diag2(:,1:5); 
[Ki]=Factor_K Long_Efectiva_Diag(coseno,nr,nudcol,X,Y,NI,NJ,L,ELEM); 
%Factor K de longitud efectiva considerando diagonales 

% Capacidad e los elementos de acero de la estructura 
[phi_Pn,phi_Vn,phi_Mn,Longitudes]=Capacidad_Elementos_Acero_Diag_Excen 
tricas(Es,Fy,Ry,Z_ELEM_COMP,Seccion2,L,Elem_Tipo_Prop2,Sec_VG_COL,Di 
mensiones_Diag,Ki,coseno,num_viga,Inic_Fin_Enlace); 

% Arriostrar a los elementos con una longitud menor al minimo entre (Lp, Lb) 





for i=1:nudcol 
D_C_Axial_omega(i,1)=max([abs(Envolvente_omega[1,1)(i,1)); 
abs(Envolvente_omega{1,2}(i,1))])/ phi_Pn(i,1); % Nudo inicial 
D_C_Axial_omega(i,2)=max([abs(Envolvente_omega{1,1}(i,4)); 
abs(Envolvente_omega{1,2}(i,4))])/ phi_Pn(i,1); % Nudo final 
end 
dibujo _Demanda_Capacidad_Axial_Columnas_Sobre_resistencia(X,Y,NI,NJ,D_C_ 
Axial_omega); 


for i=1:nuVG_COL_DIAG 

D_C_Axial(i,1)=max([abs(Envolventef{1,1}(i,1)); abs(Envolvente(1,2)(i,1))])/ 
phi_Pn(i,1); % Nudo inicial 

D_C_Axial(i,2)=max([abs(Envolvente(1,1)(1,4)); abs(Envolvente{1 ,2}(i,4))])/ 
phi_Pn(i,1); % Nudo final 


D_C_Corte(i,1)=max([abs(Envolvente{1,1}(i,2)); abs(Envolvente{1,2}(i,2))])/ 
phi_Vn(i,1); % Nudo inicial 

D_C_Corte(i,2)=max([abs(Envolvente{1 ,1}(i,5)); abs(Envolvente{1 ,2}(i,5))])/ 
phi_Vn(i,1); % Nudo final 


D_C_Flexion(i,1)=max([abs(Envolvente{1 ,1}(i,3)); abs(Envolvente{1 ,2}(i,3))])/ 
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phi_Mn(i,1); % Nudo inicial 
D_C_Flexion(i,2)=max([abs(Envolvente{1 ,1}(i,6)); abs(Envolvente(1,2)(1,6))])/ 
phi_Mn(i,1); % Nudo final 


D_C_Pn_vf(i,1)=max(D_C_Axial(i,1),D_C_Axial(i,2)); 
D_C_Vn_ví(i,1)=max(D_C_Corte(i,1),D_C_Corte(i,2)); 
D_C_Mn_ví(i,1)=max(D_C_Flexion(i,1),D_C_Flexion(i,2)); 


if D_C_Pn_vf(i,1)<0.2 
D_C_Flexo_Compr(i,1) = D_C_Pn_vf(i,1)/2+D_C_Mn_vf(i,1); 
else 
D_C_Flexo_Compr(i,1) = D_C_Pn_vf(i,1)+8/9*D_C_Mn_vf(i,1); 
end 
end 


%% Dibujos de Demanda/Capacidad 
dibujo _Demanda_Capacidad_Flexo_Compresion(X,Y,NI,NJ,D_C_Flexo_Compr); 
dibujo _Demanda_Capacidad_Cortante(X,Y,NI,NJ,D_C_Vn_vf); 





Figura 39 Código de CEINCI-LAB 


En la Figura 40 se puede observar las relaciones Demanda versus 
Capacidad de los elementos, en la Figura 40a se considera únicamente las 
acciones axiales en columnas con la demanda de la envolvente del análisis con el 
sismo amplificado por la sobre resistencia. En la Figura 40b se aprecia el efecto de 
flexo-compresión en los elementos y en la Figura 40c se considera el corte en los 
elementos, para estos análisis se considera carga muerta, viva, sismo sin 
amplificar con sus respectivas combinaciones de carga. 


Verificación de Demanda/Capacidad Axial con sobre resistencia 












D/C = 0.038 D/C = 0.023 


DIC = 0.12 DIC = 0.064 





DIC = 0.22 D/C = 0.11 





D/C = 0.35 D/C = 0.15 








[0.0 | 0.5 0.7 0.9 0.95 
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Verificación de Demanda/Capacidad *Flexo-Compresion 
DIC = 0.34 Orc = 0.35 Orc = 0.22 Ore = 0.22 OC = 0.36 DIC = 0.34 
OC = 0.21 
OC «0.15 DIC = 0.062 DIC = 0.062 DIC 013 
DIC + 0.091 DIC =0091 
, E cel yO = 022 >= 022 5 =d == 
DIC=04 DIC=04 OC = 0.22 cada pIc=022 DIC =a DIC = 04 
OC #015 orc #011 DIC*011 DIC «0.45 
DIC = 0.27 DC +027 
DIC =04t DIC=041 DIC=022 z 5 DC=022 DIC = 0.41 DIC = 41 
Cap 
DC = 0.17 DC = 0.25 DC +0235 DC = 0.17 
DIC = 0.34 DC 03 
DIC = 0.99 oc=04 oc=021 orc = 0.21 OC =04 DIC=0,39 
DC = 04 
DC=02 DIC =04 DIC =04 DC = 0.2 
DIC = 0.37 DIC=037 
0.0 0.5 0.7 0.9 0.95 
Verificación de Demanda/Capacidad "Corte 
OCc+02 Oc =02 OC +015 ~~ OC = 0.15 Orc = 0.2 Oc =0.2 
DT O19 
DIC = 008 Dro = 0.027 DY = 0.027 nC = 0.08 
bc =0 nc =o 
DIC=028 DIC =028 DIC = 0.19 oc=019 OVO = 0.28 Dro = 0.28 
OC = 0.38 
DS + 0.065 DIC = 0.094 DIC = 0.034 DIC = 0.085 
DiC = 0 pic =o 
DO + 0.26 DC #025 DIC 018 pre = 053 Oc +016 DC = 1.26 DC “0.28 
OC = 0082 DIC = 0.048 DIC = 0.048 DIC = (1092 
DIC =0 otc =0 
OC +0 OC + 0.28 OC e017 Oc +01? DC» 0.28 DIC” 028 
DC = 066 
DC = 0073 Dic = 0.062 DIC = 0.062 DIC = 0.073 
DCs oc 
0.0 0.5 0.7 0.9 0.95 


(c) 


Figura 40 a) Demanda vs Capacidad axial; b) Demanda vs Capacidad a Flexo- 


Compresión; c) Demanda vs Capacidad a Corte 


Para el diseño por capacidad se calcula el corte máximo probable con la 
ecuación 62 y se aplica una fuerza igual a 2 veces este valor en el nudo inicial del 
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elemento de enlace como se indica en el código de CEINCI-LAB de color rojo, en 
la Figura 41. Al aplicar las fuerzas de esta forma se puede simular que el elemento 
de enlace alcanza su máxima capacidad de corte. Se calculan las fuerzas en los 
elementos con estos estados de carga. 


Vur = 1.25 - Ry Van (62) 








Codigo del programa principal, parte 15: 








%% Diseño por capacidad 

% Nudo inicial y final de cada enlace 
Nudos_Enlace=nonzeros(NJ2); 

Num_inicio_ Enlaces = nuVG_COL_DIAG-num_Enlaces+1; 
Num_fin_Enlaces = nuVG_COL_DIAG; 


for i=1:num_Enlaces 
NI_enlace(i,1)=Nudos_Enlace(2*i-1); 
NJ_enlace(i,1)=Nudos_Enlace(2*i); 
Vn_Enlaces(i,1)=phi_Vn(Num_inicio_Enlaces+i-1)/0.9; 
Mn_Enlaces(i,1)=phi_Mn(Num_inicio_Enlaces+i-1)/0.9; 


V_MP(i,1) = 1.25*Ry*Vn_Enlaces(i, 1); 


M_MP(i,1) = V_MP(i,1)*L(Num_inicio_Enlaces+i-1)/2; 
Corte_NI(i,1)=-V_MP(i,1); 


Corte_Enlace(2*i-1 uve MP(i,1); 
Corte_Enlace(2*i,1)=V MEG 1); 
Momento_Enlace(2*i-1,1)=-M_M 
Momento_Enlace(2*i,1)=-M_MP(i, 

end 


%% Suponiendo sismo en Sentido de izquierda a derecha (Positivo) 
njc=num_Enlaces; 
F_Cap_P =[NI_enlace zeros(num_Enlaces,1) 2*Corte_NI 
zeros(num_Enlaces,1)]; % F = [Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=0; % Numero de miembros cargados 
Fm=0; %Elem carg, carga, código, elem a gener, incr numero elemento 
[Q_Capacidad,Q2_Capacidad]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F_ 
Cap_P,Fm,datos,nr,NI,NJ);% 
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Capacidad=KT1T\Q_Capacidad; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Numero del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 


[Fuerzas_Elementos_Capacidad_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELEM_A 
CP,L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Capacidad,Q2_Capacidad(Inic 
_Fin COL _VG(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 
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nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

ELEM_ACD2=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 

RES_ACD=0; 

% Fuerzas y momentos finales en las diagonales 


[Fuerzas_Elementos_Capacidad_DG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACDZ2, 
L_ACD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Capacidad,Q2_Capacidad(Inic_ 
Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro, RES_ACD); 


% Fuerzas en enlace 
nmiembro=1; 
RES_E=0; 


[FA_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cosE,VCE,Es,q_Capacida 
d,Q2 Capacidad(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ,nmiem 
bro, RES_B); 


% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_Capacidad_P=[Fuerzas_Elementos_Capacidad_COL_VG;Fu 
erzas_Elementos_Capacidad_DG;FA_EN]; 


%% Suponiendo sismo en Sentido de izquierda a derecha (Negativo) 
F_Cap_N = [NJ_enlace zeros(num_Enlaces,1) 2*Corte_NI 
zeros(num_Enlaces,1)]; % F = [Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=0; % Numero de miembros cargados 
Fm=0; %Elem carg, carga, código, elem a gener, incr numero elemento 
[Q_Capacidad,Q2_Capacidad]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F_ 
Cap_N,Fm,datos,nr,NI,NJ);% 
[Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_Capacidad=KT1T\Q_Capacidad; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Numero del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 


[Fuerzas_Elementos_Capacidad_COL_VG]=fuerzas_acero_arm1_V2(ngl,ELEM_A 
CP,L_ACP,seno_ACP,coseno_ACP,VC_ACP,Es,q_Capacidad,Q2_Capacidad(Inic 
_Fin_COL_VG(1,1):Inic_Fin_COL_VG(1,2),:),CG,nr,NI,NJ,nmiembro,RES); 


nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

ELEM_ACD2=ELEM_DIAG(Inic_Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:); 

RES_ACD=0; 

% Fuerzas y momentos finales en las diagonales 


[Fuerzas_Elementos_Capacidad_DG]=fuerzas_acero_arm1_V2(ngl,ELEM_ACDZ2, 
L_ACD,seno_ACD,coseno_ACD,VC_ACD2,Es,q_Capacidad,Q2_Capacidad(Inic_ 
Fin_DIAG(1,1):Inic_Fin_DIAG(1,2),:),CG,nr,NI,NJ,nmiembro,RES_ACD); 
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% Fuerzas en enlace 
nmiembro=1; 
RES_E=0; 


[FA_EN]=fuerzas_acero_arm1_V2(ngl,ELEME,LE,senE,cosE,VCE,Es,q_Capacida 
d,Q2 Capacidad(Inic_Fin_Enlace(1,1):Inic_Fin_Enlace(1,2),:),CG,nr,NI,NJ,nmiem 
bro, RES_B); 


% Fuerzas y momentos finales en toda la estructura 
Fuerzas_Elementos_Capacidad_N=[Fuerzas_Elementos_Capacidad_COL_VG;Fu 
erzas_Elementos_Capacidad_DG;FA_EN]; 


% Fuerzas y momentos finales sismic@s positivo=P y negativo=N en la estructura 
con diagonales 
Fuerzas_Elementos_Capacidad{1,1}=Fuerzas_Elementos_Capacidad_P; 
Fuerzas_Elementos_Capacidad{2,1}=Fuerzas_Elementos_Capacidad_N; 





Figura 41 Cédigo de CEINCI-LAB 


Con los estados de carga del diseño por capacidad, las cargas vivas y 
muertas, se procede a realizar las combinaciones de carga y determinar las 
envolventes de fuerzas máximas y mínimos en los elementos que corresponden a 
la Demanda (D). Anteriormente se determina la capacidad de los elementos y con 
estas variables es posible calcular la relación Demanda versus Capacidad (D/C) 
para cada elemento, el código se ilustra en la Figura 42. 








Código del programa principal, parte 15: 








tina COMBINACIONES DE CARGA CAPACIDAD „ncccsscr 
[Combinaciones_Capacidad, 
Envolvente_Capacidad]=Combinacion_Cargas(Fuerzas_Elementos_CM,Fuerzas_ 
Elementos_CV,Fuerzas_Elementos_Capacidad); 


for i=1:nuVG_COL_DIAG 
D_C_Axial_Cap(i,1)=max([abs(Envolvente_Capacidad{1,1}(i,1)); 

abs(Envolvente_Capacidad{1,2}(i,1))])/ phi_Pn(i,1); % Nudo inicial 
D_C_Axial_Cap(i,2)=max([abs(Envolvente_Capacidad{1,1}(i,4)); 

abs(Envolvente_Capacidad{1,2}(i,4))])/ phi_Pn(i,1); % Nudo final 


D_C_Corte_Cap(i,1)=max([abs(Envolvente_Capacidad{1,1}(i,2)); 
abs(Envolvente_Capacidad{1,2}(i,2))])/ phi_Vn(i,1); % Nudo inicial 

D_C_Corte_Cap(i,2)=max([abs(Envolvente_Capacidad{1 ,1}(i,5)); 
abs(Envolvente_Capacidad{1,2}(i,5))])/ phi_Vn(i,1); % Nudo final 


D_C_Flexion_Cap(i,1)=max([abs 
abs(Envolvente_Capacidad{1,2}(i,3 
D_C_Flexion_Cap(i,2)=max([abs 
abs(Envolvente_Capacidad{1,2}(i,6 


Envolvente_Capacidad{1 ,1}(i,3)); 
3)/ phi_Mn(i,1); % Nudo inicial 
Envolvente_Capacidad{1 ,1}(i,6)); 
)})/ phi_Mn(i,1); % Nudo final 


SS pr A pr, 
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D_C_Pn_vf_Cap(i,1)=max(D_C_Axial_Cap(i,1),D_C_Axial_Cap(i,2)); 
D_C_Vn_vf_Cap(i,1)=max(D_C_Corte_Cap(i,1),D_C_Corte_Cap(i,2)); 
D_C_Mn_vf_Cap(i,1)=max(D_C_Flexion_Cap(i,1),D_C_Flexion_Cap(i,2)); 


if D_C_Pn_vf_Cap(i,1)<0.2 
D_C_Flexo_Compr_Cap(i,1) = 
D_C_Pn_vf Cap(i,1)/2+D_C_Mn_vf_Cap(i,1); 
else 
D_C_Flexo_Compr_Cap(i,1) = 
D_C_Pn_vf Cap(i,1)+8/9*D_C_Mn_vf_Cap(i,1); 
end 








end 





Figura 42 Código de CEINCI-LAB para ..... 


Los gráficos de demanda capacidad para los estados de carga del diseño 
por capacidad, se presentan en la Figura 43a para flexo axial; para la relación 
demanda versus capacidad a corte, se ilustra en la 43b y el código que permite 
realizar estas gráficas se muestra en la Figura 43c. En estas Figuras se corrobora 
que las secciones de las columnas, vigas y diagonales son adecuadas con 
relaciones demanda versus capacidad menores a 1 y para los enlaces se calcula 
valores superiores a 1 debido a que se espera una mayor incursión en el rango 
inelástico. 


Verificación de Demanda/Capacidad "Flexo-Compresión 


DIC=0.27 DIC=0.3 DC=05 = qg O ne D/C = 0.3 DIC = 0.27 
DIC = 0.13 DIC = 0.13 DIC = 0.13 DIC = 0.13 
DIC = 0.84 DIC = 0.84 
DIC=0.35 DC=0.33 ii 13 DiC = 0:69 DiC = 0.33 D/C = 0.35 
DIC = 0.16 DIC = 0,37 DIC = 0.37 DIC = 0.16 
DIC = 0.8 D/C = 0.8 
DIC =0.% O/C =0.33 te da OC=0.33 DC =0.36 
DIC = 0.16 DIC = 0,53 DIC = 0.53 D/C = 0.16 
DIC = 0.81 DIC = 0.81 
DIC=0:%3 DC=0.33 DIC = 0.61 DIG=1 DIC=0:61 DIC=0.33  D/C=0.36 
D/C = 0.19 DIC = 0.73 DIC = 0.73 DIC = 0.19 
DIC = 0.76 DIC = 0.76 


Too | 0.5 0.7 0.9 0.95 
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Verificación de Demanda/Capacidad "Corte 


DIC = 0.2 BIC = 0.2 BC = 0.32 DC = 0.32 DIC =0.2 DIC = 0.2 
DIC = 1,7 
DIC = 0.08 DIC = 0,063 DIC = 0.063 DIC = 0.08 
DIC = 0.01 D/C = 0.01 
DIC=0.28  DIC=0.28 DICH O31 747 WE= 09 D/C=0.28 DIC =0.28 
DIC = 0.09 D/C = 0,088 DC = 0.088 DIC = 0.09 
DIC =0 DIC=0 
DIC = 0.28 DIC = 0.26 Beets. ja Ee 8S DC 0.28 BIC = 0.28 
VC = 
D/C = 0,082 D/C = 0,052 D/C = 0.052 DIC = 0.082 
DIC=0 DC =0 
DIC =0.28 D/C =0.28 DIC = 0.26 Dias y ¿DIC = 0.26 DC=0.28 DIC =0.28 
D/C = 0.068 D/C = 0.039 D/C = 0.039 D/C = 0.068 
DIC=0 DIC=0 
0,0 0.5 0.7 0.9 0.95 


(b) 








Codigo del programa principal, parte 15: 








%% Dibujos de Demanda/Capacidad 
dibujo_Demanda_Capacidad_Flexo_Compresion(X,Y,NI,NJ,D_C_Flexo_Compr_C 


ap); 
dibujo _Demanda_Capacidad_Cortante(X,Y,NI,NJ,D_C_Vn_vf_Cap); 








(c) 
Figura 43 a) Demanda vs Capacidad axial; b) Demanda vs Capacidad a Flexo- 
Compresión; c) Código de CEINCI-LAB 


Para finalizar el diseño por capacidad se debe calcular la relación columna 
fuerte y viga débil en base a los momentos plásticos nominales y verificar que se 
cumpla la relación de la ecuación 21. 








Código del programa principal, parte 21: 








%% Verificacion de la Conexion Viga - Columna 
[sum_M_col,sum_M_vg,Capacidad_Vg_Col]=conexion_viga_columna_V2(Fy,Fu,c 
oseno,nr,nudcol,X,Y,NI,NJ,L,Z_ELEM,Seccion,Elem_Tipo_Prop2,Sec_VG_COL,F 
uerzas_Elementos_CM,Fuerzas_Elementos_CV,Envolvente_omega); 

dibujo _Capacidad_Viga_Columna(X,Y,NI,NJ,Capacidad_Vg_Col) 





(a) 
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Verificación de Capacidad "Viga/Columna 











Vg/Col = 0.63 Vg/Col 54.1 Vg/Col = 1.1 Vq/Col = 0.63 

Wg/Col = 0.43 Vg/Col 50.82 Va/Cal = 0.82 Wg/Col = 0,43 
E 

Vg/Col = 0.44 Vg/Col 30.9 Va/Col = 0.9 Vg/Col = 0.44 
£ yq x 

Vg/Col = 0.47 Wg/Cal s4 Mg/Caol = 1 Wg/Col = 0.47 


oo] 05 07 Ao O 
(b) 


Figura 44 a) Código de CEINCI-LAB; b) Relaciones Viga. 


6. COMENTARIOS Y CONCLUSIONES 


Se expone las nuevas funciones del sistema de computación CEINCI- 
LAB para análisis y diseño de pórticos de acero con arriostramientos excéntricos. 
Se ilustra la metodología de análisis y diseño junto con la secuencia de 
programación. 


Para el prediseño de los elementos que componen al pórtico dual con 
arriostramientos excéntricos se ilustra una metodología basada en la capacidad 
del elemento de enlace, de esta forma se determina las secciones de vigas fuera 
del enlace, arriostramientos y columnas. 


Las funciones de CEINCI-LAB permiten evaluar la relación Demanda vs 
Capacidad de los elementos de un pórtico dual de manera visual, se presenta 
estas relaciones en función de colores, esto es fácil de asimilar por el usuario. 


7. RECOMENDACIONES 


Se recomienda en investigaciones futuras analizar cómo son conducidas 
las cargas gravitacionales y laterales a las cimentaciones. Es importante destacar 
que en muchos casos las cargas gravitacionales no son suficientemente grandes 
como para contrarrestar el levantamiento del pórtico excéntrico debido a las 
cargas laterales. El diseño de la cimentación es importante debido a que un 
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pórtico arriostrado excéntricamente no será capaz de alcanzar un 
comportamiento dúctil si existen levantamientos en la base. 
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